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REFERENCE 


a  Multiple-durst  Venicle  tail  Probability  Model, 
where  tne  Probability  or  Igniting 
Spillea  Flammable  Fluids  is  Time  Dependent 


i.  INTRODUCTION 


Military  vehicles  form  a  very  important  class  of  targets  against  which 
kill  probabilities  must  be  computed  for  various  weapons.  There  are  generally 
two  types  of  damage  wmcn  can  result  in  a  sill  of  a  vehicle.  These  are  (i)  a 
specified  type  or  mechanical  damage,  an d  (2)  fire.  A  fire  can  generally  be 
started  either  as  the  instantaneous  resuit  of  puncturing  some  component 
containing  flammable  fluid  with  a  nigh  energy  projectile,  or  oy  igniting  such 
fijid  wnicn  was  spiked  out  not  ignited  by  earlier  punctures.  Since  the 
amount,  and  perhaps  cne  location,  of  cne  spilled  fluid  would  generally  depend 
on  cne  numDer  and  times  of  occurrence  of  previous  punctures,  tne  probability 
of  igniting  such  spilled  fluid  might  also  depend  on  these  factors.  An  ear  lie 
report  (Reference  I)  developed  a  itiii  probability  model  including 
consideration  jf  this  dependence.  In  that  report  it  was  assumed  that  there 
was  only  or.e  .<inc  of  flammable  fluid  .n  cne  vehicle  (fuei),  ana  it  was  a!_  in 
a  single  location,  in  this  report  there  can  be  a  number  of  different  .-tines  of 
flammable  fluids,  a  number  of  different  fluid  locations,  and  up  to  two  "sinKS 
into  which  spilled  fluid  collects.  A  third  cause  of  fire  treated  in  this 
report  but  not  in  the  ear.ier  one  is  tne  detonation  or  ignition  of  stowed 
ammunition  if  it  is  struct  by  a  projectile.  The  formulas  ir.  this  report, 
then,  include  tr.ose  of  the  earlier  report  as  special  cases,  and  the  earlier 
report  can  be  considered  obsolete  and  superceded  by  this  one. 

The  model  to  oe  presented  is  not  a  venicle  vulnerability  model,  out 
merely  describes  a  method  for  putting  together  certain  basic  vulnerability- 
related  probabilities  to  obtain  the  overall  probability  tnat  the  venicle  is 
Killed.  The  basic  probabilities,  which  are  related  to  various  aspects  of 
vehicle  vulnerability,  are  assumed  to  be  available  from  standard  vehicle 
vulnerability  and  fluid  igntior.  mcceis,  and  are  not  generated  in  cne  present 


moce i . 


j.here  are  four  basic  damage  events  chat  may  occur  when  a  burst  of  rounas 
is  fired  at  a  vehicle.  These  are  (l)  the  vehicle  may  receive  disabling 
mecnamcal  damage,  (2)  some  component  containing  flammable  fluid  may  be 
punctured,  with  with  an  associated  spillage  of  that  fluid,  possibly  causing  an 
immediate  fire,  (3)  stowed  ammunition  may  be  struck,  by  a  projectile,  possible 
causing  an  immediate  fire,  and  (4)  flammable  fluid  which  was  spilled  but  not 
ignited  by  earlier  bursts  may  be  ignited  by  the  burst  under  consideration. 

Two  conditions  of  dependence  among  these  basic  events  are  treated  in  this 
report.  The  first  is  that  of  complete  independence  among  the  events,  whereby 
the  probability  of  occurrence  of  any  one  of  them  on  a  given  burst  does  not 
depend  on  whetner  any  of  tne  others  have  occurred  on  that  burst.  This  case 
mignt  be  appropriate  for  the  situation  where  the  number  of  rounds  in  the  burst 
is  large.  Tne  second  condition  is  tnat  of  complete  dependence  among  the 
events,  whereby  the  occurrence  of  any  one  of  them  precludes  the  occurrence  of 
any  or  the  others  on  the  same  burst.  This  case  might  be  appropriate  when 
consecutive  single  rounds  are  fired  at  the  vehicle,  that  is,  when  each  burst 
consists  of  only  a  single  round. 

2.  0  Lrll'v'ATiuN 

t^  be  tne  time  at  whicn  tne  rounds  from  the  i  burst  arrive  at  tne 
-arget  vehicle.  It  is  assumed  that  all  the  rounds  of  a  given  burst  arrive  at 
tne  vehicle  at  the  same  time,  but  the  time  varies  from  burst  to  burst.  This 
is  not  a  restrictive  assumption,  because  a  burst  that  extends  over  a 
significant  interval  in  time  can  be  subdivided  into  shorter  bursts,  even  to 
single-round  bursts,  if  desired,  with  each  shorter  burst  assumed  to  have  all 
its  rounds  arrive  at  tne  same  time. 

The  condition  of  the  vehicle  at  any  time  can  be  expressed  using  tne 
following  vector-like  notation: 

v  -  (F,G,g), 

where  F  describes  the  condition  of  the  vehicle  relative  to  the  puncture  of 
components  containing  flammable  fluids  which  will  drain  into  tne  first  sin*,  J 
describes  tne  condition  or  the  vehicle  rexative  to  the  puncture  of  components 


containing  fiammaoie  fluids  whicn  will  drain  into  the  second  sick,  and  g 
descriDes  the  condition  of  tne  vehicle  relative  to  disabling  mecnanical 
damage.  Specifically,  if  the  vehicle  is  not  on  fire,  F  is  tne  set  of  times  at 
which  punctures  of  the  fluid  locations  have  occurred  which  cause  fluid  to 
drain  into  tne  first  sink,  and  G  is  the  set  of  times  at  wnich  punctures  of  the 
fluid  locations  have  occurred  wnich  cause  fluid  to  drain  into  the  second  sin*. 
If  no  punctures  have  occurred,  the  symbol  0  (the  empty  set)  will  be  used  for 
either  F  or  G.  If  the  vehicle  has  suffered  disabling  mechanical  damage, 
g  ■  1;  if  nor,  g  -  0.  If  the  vehicle  is  on  fire,  the  special  vector  (*,g) 
will  ne  used.  Therefore,  using  this  notation,  a  completely  undamaged  vehicle 
would  oe  aenoted  by  tne  vector  (0,0,0). 

Berore  any  computations  can  be  made,  several  basic  probability  taoles  or 
functions  must  be  provided.  These  relate  to  the  vulnerability  of  the  vehicle, 
the  characteristics  of  the  ammunition  being  fired  at  the  vehicle,  tne  delivery 
accuracy  of  the  weapon,  and  the  number  of  rounds  in  a  burst. 


*  probability  tnat  a  puncture  of  the  j  flammable  fluid 

th 

location  occurs  on  the  i  burst. 


?  i  j ) 
io 


conditional  probability  that  tne  i  burst  causes  a  rluid 

th 

fire,  given  tnat  it  has  punctured  the  j  fluid  location. 

This  type  of  fire  will  be  referred  to  as  a  Type  I  fire,  and  is 
distinguished  from  one  caused  by  a  burst  igniting  fluid  wnich 
was  spilled  but  not  ignited  by  rounds  of  a  previous  burst. 


P1  (k) 
11 


th  th 

probability  of  perforating  the  K  sink  on  the  i  burst. 


P  (1,0 
12 


tn 

conditional  probability  that  the  i  burst  ignites  fluid  that 

was  spilled  into  sink  1  by  the  previous  punctures  identified 

th 

in  tne  set  F,  given  that  sink  1  was  perforated  on  the  i 
burst . 


P  »;  =  conditional  probability  mat  tne  i  burst  ignites  fluid  that 

id 

was  spined  into  sinx.  2  by  the  previous  punctures  identified 

tn 

in  tne  set  G,  given  tnat  sink  2  was  perforated  on  the  i 
burst. 


A  fire  resulting  from  the  ignition  of  previously-spilied 
fluids  is  called  a  Type  II  fire.  - 

Since  there  are  many  possible  compositions  of  tne  sets  F  and  G 

for  each  i,  tnere  will  be  as  many  values  of  P1  (1,F)  and 

i  ^  ^ 

P  (.2,G }.  Since  F  anc  G  are  subsets  or  the  times 
12 


i-i  i-i 

and  tnere  are  2  possible  subsets,  there  are  2  possible 

i-i 

F's  and  2  possiDxe  G's.  For  example,  if  i  «  4,  there  are 
eight  possible  F's  and  eight  possible  G's.  These  are: 


C1C2C3 


probability  tnat  the  i  burst  causes  disabling  mecnanical 
damage . 


probability  that  stowed  ammunition  is  detonated  or  ignited, 

th 

causing  a  fire,  on  tne  i  burst. 


There  are  several  other  probabilities  derived  from  those  listed  above 

that  are  required.  Ihe  first  is  the  proDability  of  puncturing  one  or  more  of 

i 

the  flammable  fluid  locations  that  drain  into  each  sinx.  If  P  (.x)  denotes 

1 

this  proDability  for  sink,  k,  then 


?  (O 
1 


(Dependent  nvencsj 


P  Cj) 

•i 


j  in  K 

k 


or 


pV) 

1 


j  in  & 


(Independent  Events) 


wnere  K  denotes  the  set  of  fluid  locations  whicn  drain  into  sink  k.  The 

k 

second  of  these  derived  probabilities  is  the  probability  of  causing  a  Type  1 

fire  in  one  or  more  of  tne  flammable  fluid  locations  tnat  dram  into  each 

sink.  if  P^k)  denotes  this  probability  for  sink  k,  then 

5 


P1(k) 

5 


(j) 


(Dependent  Events) 


j  in  K 

k 


P1(it)  =  1-/1  [1  -  P1(j)  P1  (j)) 

5  /  |  9  10 

jink 
J  k 


(Independent  Events) 


In  either  the  Dependent  Event  case  or  tne  Independent  Event  case,  the 

probability  that  fluids  are  spilled  into  sink  k  without  causing  a  Type  I  fire 
i 

denotea  P  (k),  is  given  by 
d 


P1 (k)  -  P1 (k)  -  p1 (k ) 

8  1  5 


The  final  probability  derived  from  those  shown  above  is  the  probability  of  a 

iype  II  fire,  that  is,  the  probability  that  one  or  both  sinks  are  perforated 

and  tne  fluids  spilied  into  tnose  sinks  on  tne  set  of  previous  bursts 

i 

identified  by  the  sets  F  ana  0  are  ignited.  If  P  (F,G)  denotes  this 


probability,  then 


?  (F,G)  =  P  (i)  P  (1,F)  +  p1  (2)  ?  (2,G)  (Dependent  Events) 

7  11  12  11  12 

or 

P^F.G)  ■  1  -  [  1  -  P*  (1)  P  (1,F)][1  -  (2)  (2,G)J  (Independent  Events) 

7  11  12  11  i2 

For  any  of  cnese  probabilities,  Q,  witn  the  same  subscripts,  superscripts  and 

arguments,  will  be  used  to  represent  1  -  P.  For  example, 

if  =  1  -  P1. 

4  4 

In  addition ,  note  tnat  if  F  *  0  or  G  »  0,  that  is,  if  tnere  have  oeen  no 

punctures  prior  to  time  t  ,  tnen  the  corresponding  P  *  0. 

i  12 

Now  let  S  denote  tne  set  of  possible  vectors  (vehicle  states)  after  the 
t  a  i 

i  burst.  For  i  *  0,  that  is,  before  tne  first  burst,  S  consists  of  tne 

0 

single  state  (0,0,0).  After  the  first  burst,  wnich  occurs  at  time  t  .  ten 

states  are  possible,  so  S  consists  of  ten  vectors.  These  are: 

(0,0, u)  which  occurs  if  the  first  burst  does  no  damage, 

(t  ,0,0;  wnicn  occurs  if  the  first  burst  punctures  one  or  more  f lamina  ole 

rluid  locations  wnich  drain  into  tne  first  sink,  without  causing 
a  Type  I  fire,  does  not  puncture  any  of  the  fluid  locations 
which  drain  into  the  second  sink.,  does  not  cause  disabling 
mechanical  damage,  and  does  not  cause  an  ammunition  fire, 

(0, t  0;  which  occurs  if  the  first  burst  punctures  one  or  more  flammaoie 
tluid  locations  which  drain  into  the  second  sink  without 
causing  a  Type  1  fire,  does  not  puncture  any  of  tne  fluid 
locations  which  dram  into  the  first  smic,  does  not  cause 
disaolmg  mechanical  damage,  and  does  not  cause  an  ammunition 


,  ,t . ,0) 

i  1 


wnicti  occurs  if  tne  first  burst  punctures  one  or  more  flammable 
fluid  locations  which  drain  into  the  first  sink  without  causing 
a  Type  1  fire ,  punctures  one  or  more  flammable  fluid  locations 
which  drain  into  the  second  sink  without  causing  a  Type  I  fire, 
does  not  cause  disabling  mechanical  damage,  and  does  not  cause 
an  ammunition  fire, 


(*;  0)  wnich  occurs  if  the  first  burst  does  not  cause  disabling 

mechanical  damage,  but  causes  either  a  Type  I  fuel  fire  or  an 
ammunition  fire , 

(0,0,1)  wnicn  occurs  if  the  first  burst  causes  disabling  mechanical 
damage  only, 


^,,0,1) 


which  occurs  if  tne  first  burst  punctures  one  or  more  flammable 
fluid  locations  whicn  drain  into  the  first  simc  without  causing 
a  Type  I  fire,  does  not  puncture  any  of  the  fluid  locations 
wnich  drain  into  the  second  sink,  causes  disabling  mechanical 
damage,  and  does  not  cause  an  ammunition  fire, 


, I) 


whicn  occurs  if  the  first  burst  punctures  one  or  more  flammaole 
fluid  locations  which  drain  into  the  second  sink  without 
causing  a  Type  I  fire,  does  not  puncture  any  of  the  fluid 
locations  which  drain  into  tne  first  sink,  causes  disabling 
mecnanical  damage,  and  does  not  cause  an  ammunition  fire, 


(VVI> 


wnich  occurs  if  the  first  burst  punctures  one  or  more  flammable 
fluid  locations  which  drain  into  the  first  sink  without  causing 
a  Type  I  fire,  punctures  one  or  more  flammable  fluid  locations 
which  drain  into  the  second  sink  without  causing  a  Type  I  fire, 
causes  disabling  mechanical  damage,  and  does  not  cause  an 
ammunition  fire,  and  finally, 


<*;  i) 


which  occurs  it  the  first  burst  causes  disabling  mecnanical 
damage,  and  also  causes  either  a  Type  1  fuel  fire  or  an 
ammunition  fire. 


Note  that  if  the  basic  damage  events  are  completely  dependent,  so  tnat  the 

burst  can  cause  at  most  one  of  the  events,  vectors  (t^,t^,0),  (t  ,0,0, 

(0»t  1),  (t  ,t^,l),  and  (*;  1)  cannot  occur  in  The  probabilities  for  the 

transition  from  S  to  S  are  shown  in  Table  1  for  the  case  of  Dependent 
0  i 

Events,  and  in  Table  2  for  the  case  of  Independent  Events. 


On  tne  second  burst,  similar  tables  of  transitions  from  the  vectors  in  S 
to  vectors  in  with  tne  associated  probabilities,  could  be  made.  These 


would  show  tnat  consists  ot  34  different  vectors, 


These  are : 


(0,0, g) 

(0,t  ,g) 

(trVg) 

(t2,crg) 

“lYY8’ 

(0.t  ,g> 

(0>t1t2»g^ 

“i'SY8* 

‘YbY8’ 

Ult2-tlt2-‘> 

( ^  3 , 0 ,  g ) 

(t  ,0,g) 
(t^t  ,0,g) 

(YYg> 

(tj.tj.*) 

(t  t  ,t  ,g) 

1  2  2  6 

where  g  =  0 

and  g  *  1 . 

On  the  third  burst,  similar  tables  of  transitions  from  the  vectors  in  i> 

t 

to  vectors  in  would  show  a  total  of  130  different  vectors.  These  are: 

3 


(0.0. g) 

(c  ,0,g) 

U2,0,g) 

U  .0, g) 

(c  1t3,0.g) 

U  t  ,0,g; 

2  J 

(t  t  ft  ,0,g) 


lWg> 

<t2,tr®> 

lWg) 

“iVV8’ 

“iW8’ 

‘Y^-V8’ 


10 

(t. 


4 

(t. 


U1C2 

U1C3 
(t  t 
2  3 

(t  t  t 
1  2  3 


*2.g) 

t  .g) 

t2,g) 

c2.g) 

c2.g) 

t2.g) 

c2.g) 

t2.g) 


(0 


(t2 

(t„ 

“lb 

“ibb 


t  .gl 
t-j.g) 
t-j.g) 
t3.g) 

Vg) 

t  ,,g) 

vg> 

cJ,g' 


1 


I 


Table  1 

Transition  Prooabilities  from  S  to  S 

0  1 

Dependent  Events 


Vector  in  Vector  in 


Probability 


(0,0,0) 

(id, 0,0) 

1  -  P  (1)  -  P  (2)  -  P  • 
1  1  4 

(0,t  ,0) 

1 

P  (2) 

8 

U  ,0,0) 

1 

P  (1) 

8 

(VY0> 

0 

(*,0) 

1  1 

P  (1)  +  P  (2)  +  P 
5  5 

(0,0,1) 

1 

P 

4 

(0,t  ,1) 

0 

U,0,1) 

0 

la  Die  2 

Transition  Probaoiiities  froa  S  to  S 

0  1 

Independent  Events 


Vector  in 


(0,0,0) 


(0,t  ,0) 


(t 

1 


urv0) 


(*,0) 


(0,0,1) 


(0,C  ,1) 


(t  ,0,1) 


(t  ,t  ,1) 
1  1 


(*,1) 


Probability 


Q1  ( 1 )  Q*(2)  Q'  Q 
1  1  4  o 

Q 1  ( 1 )  Pi(2)  Q1  Q1 
1  8  4  b 

1  1  11 

P  (1)  Q  (2)  Q  g 
8  1  *  o 

1  1  11 

P  (1)  ?  (2)  g  Q 
8  8  4  b 

Q1  [1  -  QJl)  Ql{2)  Q1] 
4  5  5b 

1  1  11 

Q  (1)  Q  (2)  P  Q 
1  1  4  b 

1  1  11 

P  (2)  Q  (1)  P  Q 
8  i  4  o 

1  1  1  i 

P  (1)  Q  (2)  P  g 
8  1  4  b 

1  1  11 

P  (1)  P  (2)  P  g 

8  8  4  b 

1  111 
P  [1  -  g  (1)  Q  (2)  y  j 
4  5  5 


(0,t  t 

y&) 

“’•YTa1*' 

.t,t , ,g ; 
i  - 

vw 

3>g) 

(crt2c3,g) 

*V2*g; 

U2*C1C 

yS) 

(t2’t2C3,g) 

U2-tT2t3,‘) 

(W 

3>g) 

(t3,C2t3’g) 

^ c  3 ,  t  ic  2C  3 » g ) 

(t  t  ,t  t 
12  1 

J,g) 

Ult2,t2t3,g) 

UT2'tT2t3'g) 

■SV8) 

(t  t  ,t  t 

1  3  1 

3>g) 

(tT3>t2t3’g) 

(SVtlt2t3>,) 

■wg> 

(t  t  ,t  t 

2  3  1 

3>g) 

(t2t3’t2S'g> 

<t2t3'tT2t3'g) 

•Vi'8’ 

(t  t  t  ,t  t 
12  3  1 

3>g) 

UT2t3’C2t3,S'1 

<tlV3’t.t2t3’g) 

and  (*,g),  wnere  g  *  0  and  g  -  1, 


The  pattern  snould  now  be  apparent.  in  general,  the  set  S  consists  of 
2i  2i  i 

2(2  -r  1)  vectors,  divided  into  two  sets  of  2  +1  vectors  eacn.  The  rirst 

set  has  g  *  0  and  the  second  set  has  g  *  1.  Each  of  these  sets  consists  of  a 

i  i 

single  vector  (*,g)  and  2  groups  of  2  vectors  each  corresponding  to  the 
numbers  of  possible  combinations  of  the  times 


t  ,  t  ,  t 

1  2  3 


,  •  .  t  . 

i 


before  generalizing  tne  probabilities  for  the  transition  from  set  o  to 

l-l 

set  S  ,  the  notation  to  be  used  will  be  described.  L  will  denote  the  set  oi 
i  i 

times  [t  ,  t  ,  t  ,  -  t  J.  Then  L  is  the  set  [t  ,  t  ,  t  ,  -  t  J 

1  2  3  i  i-i  1  23  i-i 

of  times  prior  to  t  .  F  and  G  will  denote  subsets  of  L  or  L  .  If  F  and  G 

l  i  - 1  i 

are  subsets  of  L  ,  then  F  t  t  and  G  +  t  will  denote  the  sets  consisting  of 
i-1  i  i 

tne  elements  of  F  and  tne  time  t  ,  or  the  elements  of  G  and  the  time  t  . 

i  i 

Thus,  for  example,  if  F  ■  [t  ,  t  ,  t  ],  then  F-t-t  «  |t  ,  t  ,  t  ,  t  |, 

124  6  1246 

Similarly,  if  G  ■  0,  then  G  +  t  ■  [t  ].  The  notation  "F  in  L  "  will  signifv 

6  6  i 

"for  all  subsets  F  of  L  Similarly,  the  notation  "G  in  L  "  will  signifv  "for 

i  i 

all  subsets  G  of  L  . “  With  this  notation  the  transition  probabilities,  and 
i 

later,  the  recursive  formulas,  can  be  very  concisely  presented.  Tables  3 

and  4  show  the  generalized  probabilities  of  transition  from  S  to  S  for  tne 

i-1  i 

Dependent-Events  case  and  the  Independent-Events  case  respectively.  (Recall 

that  P  (k,0)  *  0  for  all  i.) 
if 


A  A  A  -  "  -  * 


VA  A  > 


Table  j 

Transition  Probaoilities  from  S  to 

l-l 

Dependent  Damage  Events 

Vector  in 

Vector  in 

(For  F,  G  in  L  ) 

i-1 

S 

i 

S 

Probability 

(F, G, 0)  (F , G, 0)  1  -  P1(l)  -  P1(2)  -  P1(F,g) 

1  1  7 

(F+C  , G, 0)  P  ( U 

i  8 

(F , G-t-t  ,o)  p1(2) 

i  8 

(*,q;  p1  (.  i )  -r  p1  ( 2 )  +  p\f.g: 

5  5  / 

(F.G.l)  P* 


i  i 

P  -  P 

4  n 


Table  4 


■ 


Transition  Probabilities  from  8  to  S 

i-i  i 

independent  Damage  Events 


(For  F,  G  in  L  ) 
i-1 


Vector  in  Vector  in 


(F.G.O) 


(F,G,0) 


(F+t. ,G,0) 

1 

(F,G^t  ,Q) 
i 

(  F-t-t  ,  Cr+t  ,  0  ) 
i  i 

(*,0) 


(F.G.l) 


(F+t  , G, 1 ) 
i 

(f,G+t  ,1) 

i 

( F-t  , G+t  , 1 ) 
i  i 

(*,1) 


Probability  - 


Q  (1)  Q  U)  Q  (F,G)  q  Q 
117  4  6 

P1(l)  QX(2>  Q1(F,G)  Q1  Q1 

8  1/  h  b 

0X(1)  P1(2)  Q1 ( F , g )  q1  q1 
1  8  7  4  b 

P1(l)  P1 ( 2)  q1 ( F,G)  q*  Q1 
8  8  7  4  b 

Q1  [ 1-Qi( 1)  q1 ( 2)  Q1(F,G)  QLJ 
4  5  5  /  6 

QX(1)  Q1 ( 2)  Q1(F ,G)  P1  q1 
117  4  b 

Pi(l)  QX(2)  Q1(F,G)  P  Q1 
8  17  4  6 

q1  (l )  P1  C 2)  q1  ( F , G )  P1  q1 

18  7  4  6 

i  i  i  i  i 

P  (1)  P  (2)  q  (F,G)  P  q 
6  8  7  4  b 

P1  1 1— q1  <. l)  q1(2)  qX(F,G)  qXj 

6 


The  probabilities  in  Tables  3  and  4  will  be  used  to  develop  sets  of 

recursive  formulas  for  determining  the  probabilities  of  occurrence  of  the 

various  vectors  of  S  ,  given  tne  probabilities  for  the  vectors  of  S 

i  i-1 

First,  we  introduce  the  notation  P  (F,G,g)  to  denote  the  probability  of 

i 

occurrence  of  the  vector  (F,G,g)  in  the  set  S  .  To  illustrate  how  the 

recursive  formulas  are  developed,  consider  the  determination  of  P  (0,(9, 1). 

_ i 

Table  3  and  Table  4  show  that  the  vector  (0,0,1)  in  S  can  result  from  either 

i 

one  of  two  possible  transitions  from  vectors  in  S  .  These  are: 

i-1 

(1)  (0,0,0)  in  S  ,  with  only  disabling  mecnanical 

i~l  th 

damage  occurring  on  tne  i  burst,  or 


(2)  (0,0,1)  in  S  witn  no  puncture  or  the  fluid  system 

i  1  th 

and  no  ammunition  fire  occurring  on  the  i  burst, 


Therefore,  tne  probability  of  occurrence  of  (0,0,1)  in  S  is  the  sum  of  two 

i 

probabilities.  These  are: 

(1)  the  product  of  tne  probability  of  occurrence  of  (0,0,0) 

in  S  and  tne  probability  of  transition  from  (0,0,0) 
i-1 

in  5  to  (0,0,1)  in  S  ,  ana 
i-1  i 

(2)  the  product  of  the  probability  of  occurrence  of  (0,0,1) 

in  S  and  tne  probability  that  (0,0,1)  remains 
i“l  tn 

unchanged  by  the  i  burst. 


Thus,  for  the  case  of  independent  damage  events, 


P  (0,0,1)  -  Q/l)  QX(2)  q1  IP  (0,0,0)  P1  +  P  (0,0,1)] 
l  lib  i-1  4  i-1 


This  same  reasoning  has  been  applied  to  all  vectors  to  obtain  the  complete  set 

of  recursive  formulas  shown  in  Tables  3  and  6.  To  use  tnese  tables  to  obtain 

tne  probabilities  of  occurrence  or  tne  various  vectors  in  b  ,  ror  anv  i,  it 

i 

only  remains  to  specify  starting  values  P  (F,G,g)  for  the  possible  vectors  m 

3  .  but  b  consists  of  only  tne  single  vector  (0,0,0).  Thus  P  (0,0,0)  -  1 
0  0'  0 


Table  5 

Recursive  Formulas 

Dependent  Damage  Events 

(For  F,  G  in  L  ) 
i-1 


1.  P  (F,G,0)  -  P  ( F ,G , 0 )  [I  -  Pi(l)  -  P1(2)  -  P1(F,G)  -  P1  - J 

i  i-1  117  46 

2.  P  (F+t  , G, 0)  -  P  (F,G,0)  P1 ( 1 ) 

l  i  l-l  8 

3.  P  vF.i^rt  ,  u;  =  P  ( F  ,  G ,  0  )  PX  (  2) 

i  i  i-1  8 

4.  P  ( * , 0 )  -  P  (*,0)  y1  +  P  (F,G,0)  [ p1 (  i )  +  P1 ( 2)  +  P1 ( F , G)  +  P1 

i  i-i  4  i-i  557  b 

F , G  in  L 

i-1 

3.  P  (F,G,i>  =  P  ( F , G , 0 )  P1  +  P  (F,G, 1 )  [1  -  P1 ( 1 )  -  P1(2)  -  P^F.u) 

i  l-l  4  i-1  1  1  7 

b.  P  ( F+t  , G, 1 )  -  P  ( F , G , 1 )  P1(l) 
l  l  i-1  8 

7.  P  ( F ,  G+t  ,  1)  =  P  ( r  ,  G,  1 )  P1  (  2 ) 

l  i  i-1  8 

8.  P  (*,1)  »  P  (*,0)  P1  +  P  (*,1) 

i  i-1  4  i-i 


P  (F,G, 1 )  [ P1 ( 1 )  +  PX(2)  +  P1 ( F , G)  +  P1] 
i-1  557  6 


F,<j  in  L 


Table  6 

Recursive  Formulas 

Independent  Damage  Events 

(For  F,  G  in  L  ) 
i-1 


1.  P. (F.G.O)  -  P  (F.G.O)  Q1 ( 1 )  Q1(2)  Q1(F,G)  Q1  Q1 
i  i~l  117  46 


2.  P. (F+t  ,G,0)  -  P  (F,G,0)  P1(l)  Q1(2)  Qt(F,G)  Q1  Q1 
i  i  l-l  817  46 


-S.  P  (  F ,  G+t  ,  U  )  -  P  ( F ,  G,  0)  p1  (  2)  Q1  (  1 )  Q^F.G)  Q1  Q1 

i  i  i-1  8  1  7  46 


4.  P  (F+t  , G+t  ,0)  =  P  (F.G.O)  ?1(1)  ?1(2)  Q1(F,G)  Q1  q1 
iii  i-1  8  8  7  4  6 


5.  P  (*,U)  -  q  P  (*,0)  +  P.  (F.G.O)  [1  -  y  (1)  Q1 ( 2)  Q1(F,G)  u 
i  “  i-1  i“l  5  o  7  b 

F.G  in  L 


i-1 


6-  p .  ( F  i  G ,  1 )  =  Q\l)  y1(2)  yjF.G)  Q1  (P  (F.G.O)  P*  +  P  (F.u.l)J 

i  1  i  /  6  i-i  4  i-1 

7.  P  (F+t  ,G,1)  =*  P  U)  0  (2)  y"(F,G)  y1  [P  (F.G.O)  p"  +  P  (F.G.i'J 

ii  Hi/  6  i-1  4  i-1 

8.  P  ( F ,  G+t  ,  1 )  -  P  U)  ^  Q_(P,G)  Q1  (P  (F.G.O)  p1  +  P  (F,G,l)j 

ii  81/  b  i-1  4  i-1 

9.  P  (F+t  ,  G+t .  ,1)  *  P  (1)  P  (2)  q\f,G)  Q1  IP  (F.G.O)  p1  +  P  (F,G,l)j 

iii  8  a  7  6  i-1  4  i-1 


10.  P  (*,1)  -  P  (*,0)  P  +  P  (*,i) 
i  i-1  4  i-1 


+  [P  (F.G.O)  P^  +  P.  (F  ,G,  1)  ]  [  1  -  q\l)  Q1  (2)  Ql(F,G)  Q1 
i-i  4  i-l  5576 


F.G  in  L 


ana  P  (.F,G,g; 


0  for  all  other  (F,G,g). 


The  2(2  -t-  1)  states  whose  probabilities  of  occurrence  are  generated 

from  Tables  5  and  6  are  generally  not  of  individual  tactical  interest.  For 
example,  if  a  vehicle  is  afire  after  a  certain  number  of  bursts  have  been 
fired  at  it,  it  is  probably  of  no  tactical  interest  whether  or  not  the  vehicle 
has  also  suffered  disabling  mechanical  damage.  Table  7  shows  a  list  of  nine 
damage  categories  chat  might  be  of  tactical  interest,  the  list  of  vectors  that 
correspond  to  each  of  these  categories,  and  the  probabilities  of  occurrence  of 
each  of  these  categories  in  terms  of  the  recursively  defined  probabilities  in 
Tables  5  or  6.  -.sote  that  the  probability  of  occurrence  of  any  damage  category 
is  simply  tne  sum  of  the  probabilities  of  occurrence  of  tne  vectors  (states) 
tnat  comprise  that  category. 

Tne  damage  categories  identified  in  Table  7  can  easily  be  related  to  tne 

vehicle  kill  categories  familiar  to  army  vulnerability  analysts.  They  speak 

of  "M  or  F"  Kill,  meaning  the  loss  of  either  mobility  or  firepower,  and  "K" 

kill,  meaning  tne  catastrophic  loss  of  the  vehicle.  A  "K"  kill  occurs  only  as 

a  result  of  an  uncontrolled  fire.  An  "M  or  F"  kill  is  assessed  whenever 

certain  selected  components  are  subjected  to  specific  levels  of  mechanical 

aamage ,  and  is  also  assessed  wnenever  a  "K"  Kill  is  assessed.  Therefore  tne 

category  "Fire"  in  Table  7  can  be  equated  to  "K"  kill  as  long  as  the  Kind  of 

fire  representea  by  tne  input  values  of  P  and  P  is  the  same  kind  of 

10  12 

uncontrolled  fire  that  will  cause  a  "K"  Kill.  Similarly  if  the  mecnanical 
damage  wnose  probability  of  occurrence  is  input  as  P  is  that  wnicn  causes  an 

4 

"M  or  F"  Kill,  and  a  puncture  or  any  fluid  system  also  leads  to  sucn  a  Kill, 
then  categories  "Fluid  System  Puncture  Only"  and  "Either  Disabling  Mechanical 
Damage  or  Fire"  together  comprise  ”M  or  F"  kill. 


This  model  calls  for  many  values  of  P^(1,F)  and  P  (2,G).  For  example, 

when  computing  for  i  -  4,  there  are  three  previous  bursts,  leading  to  eight 

possible  F's  (seven  not  counting  F  =■  0)  for  which  values  of  P  (1,F)  are 

required.  The  same  number  of  values  of  P  (2,G)  would  be  required.  as  i 

1 2 

increases,  the  number  of  associated  F's  and  G's  increases  greatly.  (When 
i  =  b,  there  are  128  possible  r's  and  128  possible  G's.)  It  is  unliKeiy  that 


17 


Table  7 

Damage  Categories  and  Their  Probacilities  of  Occurrence 

After  i  Bursts 


Damage 

Corresponding  States 

Category 

(for 

all  F,G  in  L  ) 
i 

Probability  of  Occurrence 

Undamaged 

(0,0,0) 

P  (0,0,0) 

i 

Fluid  System 

(F.G.O) 

P . ( F , G , 0) 

1 

Puncture  on^v 

i  F ,  t> 

not  both  =  O] 

X 

F,G  not  both  =  0 

Disabling  Mecnanical 

Damage  Only 

(0,0,1) 

P  (0,0,1) 

i 

Fire  Only 

(*,0) 

P.(*,0) 

l 

Fluid  System  Puncture 

Disabling  Mecnanical 

Damage;  Do  tire 

1 

[  r  ,  C 

(F , u, 1 ) 

not  botn  =  0j 

P  (F,G, i ) 
i 

F,G  not  both  =  0 

Disabling  Mechanical 

Damage  ana  Fire 

(*,1) 

1 

Disabling  Mechanical 

Damage 

(F,C, 1 ) 

(*,1) 

P  C*,l)  +  P  ( F , G , 1 ) 

i  i 

F,  G  m  L 

i 


Fire 

(*,0) 

(*,1) 

P  (*,0)  +  P  (*,1) 
i  i 

Either  Disabling 

Mecnanical  Damage 

(F,G, 1) 

(*,0) 

?  (*,0)  T  ?  (*,  1)  ?  (F,v».  *. 

11  1 

or  Fire 

(*,1) 

F , u  in  - 

i 

fluid  ignition  phenomena  will  ever  be  understood  in  enough  derail  to  provide 

ail  of  these  probabilities .  however,  the  model  is  still  good.  Since  it 

provides  for  the  greatest  possible  detail  with  respect  to  the  puncture 

history,  any  less  detailed  history  is  imbedded  in  the  model.  For  example,  if 

the  only  available  data  give  the  probability  of  igniting  spilled  fluids  in 

terms  of  the  elapsed  time  since  the  first  puncture,  this  can  be  handled  by 

making  P  (1,F)  or  P  (2,G)  a  function  only  of  the  time  t  and- the  earliest 
12  >2  i 

time  in  the  set  F  or  G. 

3.  NUMERICAL  EXAMPLE 

Find  tne  probabilities  of  occurrence  of  the  damage  categories  listed  in 
Table  7  for  tne  PUSHOVER  vehicle  for  one,  two  and  three  bursts  fired  from  tne 
HL'LLY  weapon.  Maxe  tne  computations  for  botn  tne  case  where  tne  basic  damage 
events  are  assumed  to  be  dependent,  and  the  case  where  they  are  assumed  to  be 
independent.  Also  maxe  a  complete  listing  of  the  probabilities  tnat  the 
PUSHuVER  is  in  eacn  or  the  possible  damage  states  after  each  burst.  The  otne! 
input  conditions  are  as  follows: 

1.  The  POSriuV'ER  venicle  has  two  sinxs ,  the  crew  compartment 
and  tne  engine  compartment. 

2.  There  are  six  flammable  fluid  locations: 

a.  Right  Fuel  Tanx 

b.  Left  Fuel  Tank 

c.  Left  Fuel  Line 

d.  Right  Fuel  Line 

e.  Engine  Lubricant  Reservoir 

f.  Hydraulic  Cylinder  &  Lines 

3.  Locations  a,  b,  c  and  d  will  spill  fuel  into  the  crew  compartment 
if  tney  are  punctured  by  a  burst  of  fire  from  the  HULLY  weapon. 
Locations  e  and  f  will  spill  their  fluids  into  the  engine 
compartment  if  they  are  punctured. 


These  first  tnree  sets  or  required  information  vould  generally  be  availaole 
from  design  specifications  on  the  target  venicle. 


4.  The  probabilities  of  puncturing  each  of  the  fluid  locations 
on  each  of  the  bursts  is  given  by  this  table: 


Burst  1 

Burst  2 

Burst  3 

Kignt  Fuel  Tank 

.115 

.117 

.113 

Left  Fuel  Tank 

.  1U8 

.  109 

.  107 

Left  Fuel  Line 

.002 

.004 

.002 

kignt  Fuei  Line 

•  vJU  1 

.001 

.001 

t-ngine  Lubricant 

Reservoir 

.07  5 

.ObO 

.  0b  5 

Hydraulic  Cylinder  4  Lines 

.045 

.050 

.040 

These  probabilities  are 

obtained  from  one  of 

the  standard  veni 

vulnerability  models. 

5.  The  probabilities  or 

an  immediate  catastrophic  fire 

,  given 

that  each  of  these  locations  are 

punctured 

on  each 

of  the 

bursts  is  given  by  tnis  table: 

Burst  1 

Burst  2 

Burst  3 

Right  Fuei  Tank 

.750 

.760 

.740 

Left  Fuel  Tank 

.C>50 

.  640 

.  ob 0 

Left  Fuel  Line 

.500 

.400 

.600 

Right  Fuel  Line 

.200 

.350 

.250 

Engine  Lubricant 

Reservoir 

.050 

.060 

.070 

Hydraulic  Cylinder  i»  Lines 

.040 

.030 

.050 

These  probabilities  are  obtained  either  from  test  data  on  fluid  ignition 
or  from  a  theoretical  model  or  such  ignition. 


P. 


Tne  probabilities  that  each  of  the 
or.  eacn  of  tne  Bursts  are  given  in 


sinks 

this 


are  pe 
taole  : 


rf orated 


Crew  Compartment 
Engine  Compartment 


Burst  1  Burst  2  Burst  3 


.450  .500  .400 

.300  .350  .325 

These  probabilities  are  obtained  from  one  of  the  standard  vehicle 
vulnerability  models. 

7.  The  probaDilities  of  igniting  fluid  spilled  into  each  sink 
on  previous  bursts  given  tnat  burst  2  perforates  the  sink 
are  given  in  tnis  table: 

Fluids  Spilled  on 
Burst  i 

Crew  Compartment  .250 

engine  Compartment  .075 

Similar  probabilities  for  burst  3  are  given  in  this  table: 

Fluids  Spilled  on 
Burst  1  Burst  2  Botn  1  &  2 
Crew  Compartment  .275  .2b5  .235 

Engine  Compartment  .080  .070  .065 

These  prooacilities  are  obtained  either  from  test  data  on  fluid  ignition, 
or  from  a  tneoretical  model  of  sucn  ignition. 

8.  The  probabilities  of  disabling  mechanical  damage  on  each  of 
the  bursts  are  given  in  this  table: 

Burst  1  Burst  2  Burst  3 

Disabling  Mechanical  Damage  .150  .170  .160 
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Tnese  probabilities  are  obtained  trom  one  ot  tne  standard  venic^e 
vulnerability  models. 

9.  The  probabilities  of  igniting  stowed  ammunition  on  each  of 
the  bursts  are  given  in  tnis  table: 

Burst  1  Burst  2  Burst  3 

Igniting  Stowed  Ammunition  .100  .110  .090 

Tnese  prooaoilities  are  obtained  from  one  of  the  standard  vehicle 
vulnerability  models. 

The  input  set  corresponding  to  these  conditions  and  prepared  according  to 
the  instructions  given  in  the  comments  in  the  program  itself  (aPPENDIX  A)  is 
shown  in  Table  8. 

The  output  of  this  run  consists  of  five  pages  of  line-printer  output. 

The  first  page  (Table  9)  consists  only  of  a  recap  of  the  input  data,  labeled 
so  tnat  it  will  be  easy  tc  cnecK  it  for  correctness.  The  second  page 
(Table  10)  lists  the  probabilities  that  the  target  vehicle  is  in  each  of  the 
damage  categories  snown  in  Table  7,  after  tne  first  burst,  after  the  second 
burst,  and  after  the  third  burst.  Also  shown  are  the  more  common  Kill 
probabilities  in  terms  of  M,  F  and  X  Kills,  and  how  they  are  related  to  the 
nine  damage  categories.  Results  are  shown  for  both  the  assumption  of 
dependent  damage  events  and  the  assumption  of  independent  damage  events.  If 
we  had  asRed  for  only  one  of  these  assumptions,  this  page  would  have  had  the 
results  for  that  assumption  only.  Tables  11,  12  and  13  show  the  probabilities 
that  the  vehicle  is  in  each  of  the  many  damage  states  after  each  burst.  These 
results  may  be  omitted  by  the  choice  of  the  last  input  number. 

Eacn  page  of  machine  output  is  labeled  at  the  top  with  the  first  two 
input  lines  that  allow  for  arbitrary  user-supplied  case  descriptors.  Each 
page  of  output  is  also  labeled  witn  the  date  and  time  that  the  macnine  run  was 
made.  This  provides  a  furtner  means  of  identifying  runs  and  of  matching  input 
secs  (output  page  1)  with  output  sets  (output  pages  2-3). 
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Table  8 

Input  File  for  Numerical  Example 


Target  Vehicle  -  PUSHOVER 
Tnreat  Weapon  -  BULLY 


Taole  9 

First  Output  Page  -  Numerical  Example 


Target  Vehicle  -  PUSHOVER 
Threat  Weapon  -  BULLY 

INPUT  LIST 

Number  of  Bursts  ■  3  Number  of  Fluid  Locations  *  6  Number  of  Sinks  *  2 

Fluid  Locations  Draining  into  Sin*  1 :  1234 

Sin*  2:  5  6 


BURST 


1 

2 

3 

Prob  of  a  Fluid  System  Puncture  (P9)  Loc  1 

0.1150 

0.1170 

0.1130 

Loc  2 

0.1080 

0. 1090 

0. 1070 

Loc  3 

0.0030 

0.0040 

0.0020 

Loc  4 

0.0010 

0.0010 

0.0010 

Loc  3 

0.0750 

0.0800 

0.0650 

Loc  t> 

0.0450 

0.0500 

0.0400 

Prob  of  Type  I  Fire/ Puncture  (P10)  Loc  1 

0.7500 

0.7600 

0.7400 

Loc  2 

0.6500 

0. 6400 

0. 6o0o 

Loc  3 

0.5000 

0.4000 

O.oOOO 

Loc  4 

0. 30uu 

0. 3500 

0. 2500 

Loc  3 

0.0500 

0. ObOO 

0.0700 

Loc  6 

0.0*00 

0. 0300 

0.050'J 

Proo  of  Disabling  Mecnanical  Damage  (P4) 

0. 1500 

0. 1700 

0. 1600 

Prob  of  r ire  Due  to  ttit  on  Ammunition  (P6) 

0. 1000 

0.1100 

0. C900 

Probability  of  Perforating  SinKs  (Pli)  Sink  1 

0.4500 

0.5000 

0.4000 

Sink  2 

0.3000 

0.3500 

0.3250 

Probability  of  Igniting  Fluid  Spilled  on  Set  F 
Given  a  Perforation  of  the  Appropriate  Sink 

of  Previous 
[ P12( F )  ] 

Bursts , 

F  *  1 1 

F  -  t2  F 

-  tlt2 

Burst  2,  Sink  1  0.2500 

___ 

Burst  3,  Sink  1  0.275U 

0.2650  0 

.2850 

Burst  2,  Sink  2  0.0750 

— 

Burst  3,  Sink  2  U.08GU 

^ i  CO  o 

.  0t>5C 

This  run  was  made  Sun  Jan  12  14:2o;u_  UST  198b 
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Table  10 

Second  Output  Page  -  Numerical  Example 


Target  Vehicle  -  PUSHOVER 
Threat  Weapon  -  BLILEY 

OUTPUT  LIST 


Calculations  Made  Assuming  Dependent  Events 

NUMBER  OF  BURSTS 


1  2  3 


1. 

Undamaged 

0.40300 

0.14468 

0.06105 

t—  « 

Puncture  Only 

0.1832O 

'O 

3' 

r^. 

o 

• — 4 

• 

o 

0. 11274 

3, 

Mechanical  Only 

0. 1 50OO 

0. 1478b 

0.  1092'J 

4. 

Fire  Only 

o. 2b38Q 

0. 39289 

0.  -*1882 

5. 

Puncture,  Mechanical,  No  Fire 

— 

0.06025 

0.09310 

t>. 

Mechanical,  Fire 

— 

0. 08b39 

O.2u5o7 

7 . 

Mechanical 

0.15000 

0.29450 

0.40738 

3. 

Fire 

0.26380 

0.47927 

0.62390 

9. 

Eitner  Mechanical  or  Fire 

0.41380 

0.68739 

0.82620 

Neither  M,  nor  F,  nor  K  (1) 

0.40300 

0. 14468 

0.06105 

M 

or  F,  but  not  K  (2+3+5) 

0.33320 

0. 37b05 

0.31505 

M 

or  F  (including  K)  (2+9) 

0.59700 

0.85532 

0.93895 

K 

(8) 

0.26380 

0.47927 

0.O2390 

Calculations  Made  Assuming  Independent  Events 

NUMBER  OF  BURSTS 


1 

2 

j 

1. 

Undamaged 

0.53134 

0.26855 

0.  14551 

2. 

Puncture  Only 

0.11384 

0.12724 

0. 10506 

3. 

Mechanical  Only 

0.09377 

0.11210 

0.10003 

4. 

Fire  Only 

0.20482 

0.30971 

0.34203 

5. 

Puncture,  Mechanical,  No  Fire 

0.02009 

0.05312 

0.07223 

6. 

Mechanical,  Fire 

0.03614 

0.  12928 

0.23512 

7. 

Mechanical 

0.15000 

0.29450 

0.40738 

8. 

Fire 

0.24096 

0.43900 

0.57715 

9. 

Either  Mechanical  or  Fire 

0.35482 

0.60421 

0.74941 

Neither  M,  nor  F,  nor  K  (1) 

0.53134 

0. 2bb55 

0. 14551 

M 

or  F,  but  not  K  (2+3+5) 

0.22769 

0.29246 

0.27734 

M 

or  F  '.including  K.)  (2  +  9) 

0. 46866 

0.  7  3145 

0.85*449 

R 

(8) 

0.24096 

0.43900 

0.57715 

This  run  was  made  Sun  Jan  1*.  14: 2b:  12  EST  1986 
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Table  11 

Third  Output  Page  -  Numerical  Example 


Target  Venicie  -  PUSHOVER 
Threat  Weapon  -  BULLY 


Individual  State  Probabilities  -  Dependent  Events 


-  Burst 

1 - 

-  Burst  2  - 

-  Burst  3  - 

F 

G 

g  -  0 

g  "  1 

g  -  0 

g  “  1 

g  *  o 

g  *  1 

phi 

phi 

0. 10300  0 

.15000 

0. 14468 

0. 14786 

0.06105 

0.10920 

pni 

ti 

0.11445 

— 

0.03808 

0.01946 

0.01508 

0.01691 

phi 

t2 

— 

— 

0.04985 

0.01856 

0.0199U 

0. 018j5 

phi 

t3 

— 

— 

— 

— 

0.01424 

0. 01456 

phi 

t  It  2 

— 

— 

0.01416 

— 

0. 005o8 

0.00227 

phi 

t  It  3 

— 

— 

— 

— 

0.00375 

0.00192 

phi 

t2t  3 

— 

— 

— 

— 

0.00491 

0.00183 

pni 

t  lt2t  3 

— 

— 

— 

— 

0.0U139 

— 

tl 

pni 

0.06875 

— 

0.01609 

0.01169 

0.00502 

0.00809 

tl 

1 2 

— 

— 

0.00850 

— 

0.00246 

0.00136 

tl 

t3 

— 

— 

— 

— 

0.00158 

0.00115 

tl 

t2t3 

— 

— 

— 

— 

0.00084 

— 

1 2 

phi 

— 

— 

0.02836 

0.01056 

0.00896 

0.00956 

t2 

tl 

— 

— 

0.00805 

— 

0.0023** 

0.00129 

1 2 

1 3 

— 

— 

— 

— 

0.00279 

0.00104 

1 2 

1 1 1 3 

— 

— 

— 

— 

0. 00u7  9 

— 

t  J 

pni 

— 

— 

— 

— 

0.00974 

0.00995 

t3 

tl 

— 

— 

— 

— 

0.0025o 

0.00131 

1 3 

1 2 

— 

— 

— 

— 

0.00336 

0.00125 

1 3 

t  It  2 

— 

— 

— 

— 

0.00095 

— 

t  It  2 

pni 

— 

— 

0.00484 

— 

0.00149 

0.00077 

t  it2 

1 3 

— 

— 

— 

— 

0.00048 

— 

t  it  3 

phi 

— 

— 

— 

— 

0.00108 

0.00079 

t  lt3 

t2 

— 

— 

— 

— 

0.00057 

— 

1 2t  3 

pni 

— 

— 

— 

— 

0.00191 

0.0007 1 

1 2t  3 

tl 

— 

— 

— 

— 

0.00054 

— 

t lt2t  3 

phi 

— 

— 

— 

— 

0.00033 

— 

* 

This  run  was  made 

0.26380 

Sun  Jan  12 

14:26:42 

0.39289 

EST  1986 

0.08639 

0.41882 

0. 205u7 

Table  12 

Fourth  Output  Page  -  Numerical  Example 

Target  Venicle  -  PlSriOVt.it 
Threat  Weapon  -  al'LLY 


Individual  State  Probabilities  -  Independent  Events 


-  Burst  1  - 

-  Burst  2  - 

-  Burst  3  - 

F 

G 

g  *  o 

g  ”  1 

g  *  o 

g  “  1 

g  -  o 

g  “  1 

phi 

phi 

0.53134 

0.09377 

0.26855 

0.11210 

0.14551 

0.10003 

phi 

tl 

0.06681 

0.01179 

0.03288 

0.01373 

0.01735 

0.01193 

pni 

t2 

— 

— 

0.03678 

0.01525 

0.01948 

0.01339 

phi 

t3 

— 

— 

— 

— 

0.01554 

0.01068 

pni 

t  It 2 

— 

— 

0. 00450 

0.00188 

0.00^39 

0.00164 

pni 

t  lt3 

— 

— 

— 

— 

0.00185 

0.00127 

phi 

t2t3 

— 

— 

— 

— 

0. 00208 

0.00143 

pni 

t  It  2t  3 

— 

— 

— 

— 

0.00026 

0.00018 

tl 

phi 

0.04177 

0.00737 

0.01847 

0.00771 

0.00891 

0.00612 

t  i 

tl 

0. 00525 

0.00093 

0.00226 

0.00094 

0. OU106 

0. 0007  i 

1 1 

t2 

— 

— 

0.00253 

0.00106 

0.00119 

0.00082 

tl 

1 3 

— 

— 

— 

— 

0.00095 

0.00065 

1 1 

t  lt2 

— 

— 

0.00031 

0.00013 

0.00015 

0.000 10 

t  i 

t  it  3 

— 

— 

— 

— 

0.00011 

0.00008 

tl 

t2t3 

— 

— 

— 

— 

0.00013 

0.00009 

1 1 

t  it  2t  3 

— 

— 

— 

— 

0.00002 

0.00001 

1 2 

pni 

— 

— 

0.02163 

0.00903 

0.01048 

0.00720 

1 2 

tl 

— 

— 

0.00265 

0.00111 

0. 00125 

0.00086 

1 2 

L.  *1 

— 

— 

0.00296 

0.00124 

0.00140 

0.00096 

1 2 

t  3 

— 

— 

— 

— 

Q.Oul  i2 

0 . ObC  7  / 

1 2 

tit*. 

— 

— 

0.00036 

0.00015 

0.00017 

0.00012 

1 2 

t  it  3 

— 

— 

— 

— 

0.00013 

0. 000o9 

1 2 

t2t3 

— 

— 

— 

— 

0.00015 

0.00010 

1 2 

t  It2t3 

— 

— 

— 

— 

0. 00uo2 

0.00001 

t3 

pni 

— 

— 

— 

— 

0.01119 

0.0o7t>9 

c  3 

t  i 

— 

— 

— 

— 

0. 00133 

0. 0009 2 

t3 

1 2 

— 

— 

— 

— 

0.00150 

0.00103 

1 3 

1 3 

— 

— 

— 

— 

0.00119 

0.  uOOs*. 

t3 

t  it  2 

— 

— 

— 

— 

0.00018 

0.00013 

1 3 

tlt3 

— 

— 

— 

— 

0. OOOl 4 

0.00010 

t3 

t2t3 

— 

— 

— 

— 

0.00016 

0.00011 

t  3 

tlt2t3 

— 

— 

— 

— 

O.OOOOi 

0.00001 

t  lt2 

phi 

— 

— 

0.00149 

0.00062 

0.00071 

0.00049 

t  lt2 

tl 

— 

— 

0.00018 

0.00008 

0.00009 

0.00UU6 

t  it2 

t2 

— 

— 

0.00020 

0.00009 

0.00010 

0.00007 

t  It  2 

t3 

— 

— 

— 

— 

0.00008 

u.  uOOo  :> 

t  lt2 

t  lt2 

— 

— 

0.00002 

0.00001 

0.000C1 

0.00001 

t  lt2 

1 1 1  3 

— 

— 

— 

— 

0.00001 

O.OOOl  . 

t  it  2 

t2t3 

— 

— 

— 

— 

0. 00001 

0.  0o"J,  1 

t  lt2 

tlt2t3 

— 

— 

— 

— 

o. OOUuc 

o. Ouuol 

Table  Continued  on  Next  Page 

This  run  was  made  Sun  Jan  12  11:26:42  EST  1986 
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Fifth 

.  Output  Page  - 

Numerical  Example 

Target  Vehicle 

-  PUSHOVER 

Threat  Weapon 

-  HOLLY 

Individual  ! 

State  Probabilities  -  Independent 

Events 

(continued) 

-  LUlSt  1  - 

-  Burst  2  - 

-  Burst  3  - 

F  0 

g  *  0 

g  ’  1 

g  *  0  g  »  1 

g  -  0 

g  -  1 

tlt3  phi 

— 

— 

—  — 

0.00068 

0.00047 

tit!  t 1 

— 

— 

— 

0.00008 

0.00006 

1 1 1 3  t2 

— 

— 

—  — 

0.0U009 

0. OOuuo 

tit!  t3 

— 

— 

— 

0.00007 

0.00005 

t  It3  t It  2 

— 

— 

—  — 

0.00001 

O.OOOOI 

t  it  3  t  it  3 

“ 

— 

—  — 

0.000U1 

0.00001 

t  it3  t2t  3 

— 

— 

—  — 

0.00001 

o.ooooi 

t  It 3  t  It2t3 

— 

— 

—  — 

0.00000 

0.00000 

1 2t  3  pni 

— 

— 

—  — 

0.00081 

0.00055 

t2t3  tl 

— 

— 

— 

0.00010 

0.00007 

t2t3  t2 

— 

— 

—  — 

0.00011 

0.00007 

t2t3  t3 

— 

— 

—  — 

0.00009 

0.00006 

t2t3  t lt2 

— 

— 

—  — 

0.00001 

0.00001 

t2t3  t it3 

— 

— 

— 

0.00001 

0.00001 

t2t3  t2t3 

— 

— 

—  — 

o.ooooi 

O.OOUOl 

t2t 3  tlt2t3 

— 

— 

— 

0.0U000 

O.OOuOO 

tlt2tj  pni 

— - 

— 

—  — 

0.00005 

0.00004 

t  it  2t  3  cl 

— 

— 

— 

O.OOOOi 

0.00000 

tit2t3  t2 

— 

— 

—  — 

0.00001 

0.00001 

tlt2t3  t3 

— 

— 

— 

0.00001 

0.00000 

t  It  2t  3  t  It  2 

— 

— 

— 

0. 00000 

0.  OOUUO 

Clt2t3  t it 3 

— 

— 

— 

0.00000 

0. U00G0 

t It2t3  1 2t 3 

— 

— 

—  — 

U.OOUUO 

0.  OlOOO 

t It2t3  t  It2t3 

— 

— 

— 

0.00000 

0. 00000 

* 

0.20482 

0. 03b  14 

0. 30971  0.12928 

0.34203 

0.23512 

This  run  was  made  Sun  Jan 

12  14:26:42  EST 

1986 

The  damage  categories  whose  probabilities  of  occurrence  are  given  in 
Table  iO  fall  into  three  classes.  These  are: 

Class  1  -  Categories  that  cannot  be  entered  from  other 

categories,  only  left  in  moving  to  other  categories. 

Class  2  -  Categories  that  can  be  both  entered  from  otheTT" 

categories  and  left  in  moving  to  other  categories. 

Class  3  -  Categories  that  can  only  be  entered  from  other 
categories,  but  never  left. 

uf  tne  categories  listed,  the  only  Class  1  category  is  "Undamaged."  This  is 
tne  initial  state  of  the  target  before  any  firing.  Class  2  categories  are 
those  representing  only  partial  damage.  These  are  "Puncture  Only," 

"Mecnanical  only,"  "Fire  only,"  and  "Puncture,  Mecnanical,  No  Fire."  Class  3 
categories  are  those  representing  either  partial  or  complete  damage.  These 
are  "Mecnanical,  Fire,"  "Mechanical,"  "Fire,"  and  "Either  Mechanical  or  Fire." 
As  the  number  of  bursts  increases,  the  probabilities  associated  witn 
categories  in  these  classes  will  cnange  in  different  ways.  Probabilities  for 
Class  i  categories  will  always  decrease,  for  subsequent  bursts  will  always 
move  tne  target  out  of  the  category,  never  into  it.  Probabilities  for  Class  2 
categories  will  generally  increase  at  first,  but  finally  decrease,  and,  in  tne 
limit,  become  zero.  This  is  because  Class  2  categories  can  be  entered  irom 
Class  1  categories  or  other  Class  2  categories,  and  as  the  probabilities 
associated  with  tnese  "feeder"  classes  decrease,  there  is  less  chance  of  the 
target  entering  the  Class  2  category;  it  can  only  leave  it.  Probabilities  for 
Class  3  categories  always  increase,  approaching  unity  in  the  limit.  If  an 
infinite  number  of  bursts  are  fired  at  the  target,  it  will  with  complete 
certainty  suffer  botn  "Fire"  and  "Mechanical"  damage. 

These  trends  are  noted  in  Table  10  although  the  probabilities  for  some 
of  the  Class  2  categories  have  not  yet  started  to  decrease  after  three  bursts . 
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This  program  implements  the  metnodology  in  tnis  report  for  up  to 
tnree  oursts.  The  code  is  written  in  FORTRAN  77.  Instructions 
for  preparing  the  inputs  are  included  as  comments  preceding  each 
of  tne  "read”  statements;  all  input  is  read  in  at  the  beginning 
of  the  program.  All  numerical  input  is  read  using  " Lis t-Uirected 
RiAiT  statements.  (See  the  example  in  tne  report  text.; 

implicit  real*8  (a-h,  o-z) 
character*80  ident 1 ,  ident2 
character*60  today 

common  pl(3,2),  p3(2,2:3,3),  p4(3),  p5(3,2),  p(6,9),  p6t_3_), 

*  p9( 3,10),  P10( 3,10),  pll( 3,2) ,  p 12( 2, 2: 3,3), 

*  q i ( 3, 2 ) ,  q3( 2, 2: 3, 3) ,  q4(3),  q5(3,2),  q6(3), 

*  s(65,b),  t(65,6),  nis(2),  isinx(2,10),  pbrl(3,4) 


The  first  information  read  is  tne  date  and  time  of  tne  macnine 
run.  At  AMsAA,  where  this  program  is  implemented  on  a  VAX  11/78C 
wim  tne  L’MX  operating  system,  tnis  information  is  adaea  to  the 
input  rile  automatically  by  an  executable  script  that  inaentifies 
for  the  program  tne  input  file  name,  and  issues  the  command  to 
run.  Therefore  it  should  not  be  typed  into  the  input  file  by  the 
user.  If  tnis  program  is  implemented  by  other  users  at  other 
agencies,  this  information  snould  be  provided  in  whatever  way  is 
appropriate,  or  else  tne  "read”  statement  deleted  from  tne  code. 

read (*,8000)  today 

format (abO) 


The  first  user-supplied  input  consists  of  2  lines 
(up  to  80  columns  each)  of  arbitrary  information  to 
identify  or  explain  the  case  being  run.  Both  lines 
must  be  provided,  even  if  one  is  empty. 

read(*,8UUl)  ident  1 
reaa(*,8001)  ident! 
rormat (aBO) 


The  second  input  is  a  single  integer  telling  the  number 
of  bursts  in  the  case  being  run. 

The  maximum  allowable  number  is  three. 


read(*,*)nr 


The  tnirn  input  is  a  single  integer  telling  tne  number 
of  sinks  into  wnich  spilled  fluid  can  collect. 

The  maximum  allowable  number  is  two. 


reaa(*,*)nsinxs 


The  fourth  input  is  a  single  integer  telling  the 
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c  number  of  fluid  locations  on  the  vehicle, 

c  The  maximum  allowable  numoer  is  lU. 

c 

read(*,*)  nil 


c 

c  The  fifth  input  consists  of  nr  lines,  one  for  each  burst, 

c  The  ith  line  has  nfl  numbers,  the  jth  of  which 

c  gives  the  probabilities  p9  for  tne  jth  fluid  location, 

c  (p9  is  the  probability  of  puncturing  the  jth  fluid  location 

c  1  on  the  ith  burst.)  -  — 

c 

do  100  i  »  1,  nr 
read(*,*)(p9(i , j ) ,j*l,nfl) 
luO  continue 

c  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  — 

c  The  sixth  input  consists  of  nr  lines,  one  for  each  burst, 

c  The  ith  line  nas  nfl  numbers,  the  jth  of  which 

c  ,  gives  tne  probabilities  piO  for  the  jth  fluid  location, 

c  (plO  is  the  probability  of  type  1  fire,  given  a  puncture) 

c 

do  110  i  =  1,  nr 

read{ * , *) (p 10(i , j ),j=l,nfl) 

110  continue 

c 

c  The  seventh  input  consist  of  two  numbers;  the  first  is 

c  the  number  of  fluid  locations  which,  if  punctured,  will 

c  drain  into  tne  first  sink,  and  tne  second  number  is  the 

c  number  of  fluid  locations  wnicn,  if  punctured,  will 

c  drain  into  the  second  sinx. 

c  If  there  is  only  one  sink,  this  input  should  not  be 

c  provided:  the  program  will  automatically  assign  sinx  one 

c  to  all  locations, 

c 

it C  ns  inns  .eq.  2  )  read(*,*)(  nis(i),  i  =  1,2  ) 


c 

c  The  eighth  input  consists  of  a  mapping,  teliing  whicn 

c  fluid  locations,  if  punctured,  will  drain  into  each  sink, 

c  If  there  is  only  one  sink,  nothing  at  all 

c  needs  to  be  input.  The  program  assumes  that  all  of  the 

c  fluid  locations  drain  into  that  sink.  If  there  are 

c  two  sinxs,  you  must  input  2  lines  here;  the  first  line 

c  lists,  as  integers,  the  numbers  of  the  fluid  locations  that 

c  drain  into  the  first  sink,  and  the  second  line  lists  the 

c  numbers  of  the  fluid  locations  that  drain  into  the  second 

c  sink, 

c 

i£(  nsinxs  .eq.  2  )tnen 
do  120  i»i , 2 

read  (*,*  His  ink  (i  ,  j  ) , j*i,nis(i ) ) 

120  continue 

else 
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nis  (  1 )  =*  nr  I 
do  1  JO  j  =  1 , n 
isintcC  l,j  )=j 
continue 
endif 


Tne  nintn  input  consists  of  nsinks  lines.  The  ith 
of  these  lines  gives  pll's,  the  probabilities  that  the  ith 
sink  is  perforated  on  each  of  the  bursts. 

do  140  j  =*  1 ,  nsinks 
read(*,*)(  pll(i,j),  i-1,  nr  ) 
continue 


The  tenth  input  is  a  set  or  lines  giving  tne  proDabixities 
of  igniting  fluids  spilled  on  previous  bursts 
(pl2),  for  tne  second  and  tnird  burst. 

The  first  of  tnese  lines  gives  tne  probaDility  pl2(l,2,F) 
for  F  =  1.  (for  the  second  burst  and  tne  first  sink.) 

The  second  of  tnese  lines  gives  the  probabilities  pl2(l,3,F) 
for  F  =  1,  F  =  2,  and  F  =  1,2.  (for  the  tnird 

burst  and  the  first  sins..) 

The  third  of  tnese  lines  gives  the  probability  pi 2(2, 2, F) 
for  F  =  1.  (for  the  second  burst  and  tne  second  sinK.  ) 

The  fourth  of  these  lines  gives  the  probabilities  pl2(2,3,F) 
for  F  =  1,  F  =  2,  and  F  =  1,2.  (for  the  third 

burst  and  tne  second  sinkt.) 


if(  nr  .gt.  1  )  then 

do  160  i  *  2, nr 
m=2**(i-l )-l 
do  150  j  =  1 ,m 
P3( 2, i , j )  =*  0. 
q  3(  2 , i , j  )  *  1. 
continue 
continue 
do  190  k»l, nsinks 
do  180  i  =  2, nr 
m  *  2**(i-l)  -  1 
read(*,*)(pl2(k,i,j ),j-l,m) 
do  170  j  =  1,  m 

p3(*t ,  i ,  j  )  «  pil(i.k)  *  pl2(k,i  ,  j  ) 
q3(k,i,j )  *  1.  -  p3(k,i,j ) 
continue 
continue 
continue 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


The  eleventh  input  line  gives  the  probabilities 
or  disabling  mechanical  damage  (p4)  for 
the  first,  second,  and  third  burst. 

read(*,*)(p4(i) ,i*l,nr) 


The  twelvth  input  line  gives  the  probabilities 
of  fire  due  to  a  hit  on  stowed  ammunition  (po)  for 
Che  first,  second,  and  third  burst. 

read(*,*)(p6(i) ,i»i,nr ) 
do  200  i  =  1 ,  nr 
q4(i)  -  1.  -  p4(i) 
q6(i)  -  1.  -  po(i ) 

200  continue 


The  thirteenth  input  line  is  an  indicator  telling  wnetner 
to  calculate  using  tne  formulas  for  dependent  events  or 
tnose  for  independent  events,  input  0  for  dependent  events 
only;  input  1  for  independent  events  only;  and  input  2 
for  botn. 

reacK*,*;icor 


The  last  input  line  is  an  indicator  telling  whether  to 
print  tne  probabilities  of  occurrence  of  the  individual 
vectors  (states).  Input  1  to  print  out  these  probabilities, 
or  0  to  not  print  them. 


read(*,*)ipv 

c----------- 


c  write  out  the  entire  input  set  for  this  run. 

c 

c-----------  -  -  --  --  --  --  --  --  --  --  -- 

write (*, 8 123)  identl 
8128  format (' 1 ’adO) 

write (*, 8Gu2 )  identl 
8002  format('  'a 80,//; 

write(*,8003) 

80u3  format (31x, ' I  NPUT  LIST',) 

wrice(*,8004)  nr,  nfl,  nsinks 
8004  format(/'  Number  of  Bursts  *',i2,5x, 

*  'Number  of  Fluid  Locations  *'ii,4x, 

*  'Number  of  binks  *',i2/) 

if(nsinks  .eq.  2)then 
do  2lu  i  *  i,  nsinics 

if(  i  .eq.  I  )write ( * , 8005 )  i,  (isinK(i.j),  j  *  i,nis(i)) 
if (  i  .eq.  2  )wri te ( * , 800b)  i,  (isimt^i , j ) ,  j  «  1, nisei)) 
2 ivJ  continue 
endif 


8005  format (’  Fluid  Locations  Draining  Into  Sink  ' ,i 1 , ' : ' , i0i3) 
SOOo  format('  Sink  * ,i 1 , ' : ' , iOii) 


write ( * , 8007 ) 

8007  format (bix, '8  U  R  S  T ' , // 53x , ' 1 ' , lOx , ’ 2 ' , lOx , ' 3 ' ) 


write (*, 8008) (p9(i , 1) ,i«l, nr) 

8008  format(/'  Prob  of  a  Fluid  System  Puncture  (P9)  Loc  1', 
*  2x,3(f 11.4)) 

if (  nfl  .gt.  l)then 
do  220  j  ■  2,  nfl 
write e*,80C9)j , ep9e i , j ) ,i-l .nr ) 

8009  format e 38x , '  Loc ' ,i2, 2x, 3ef 1 1 . 4) ) 

220  continue 

endif 


write ( * ,8010)epl0ei,l),i“l,nr) 

SulO  formate/'  Proo  of  Type  I  Fire/ Puncture  (P10)  Loc  1', 
*  2x,3(fli.4)) 

ife  nfl  .gt.  l)tben 
do  230  j  *  2,  nfl 
wr i tee*, 8009 )j , eplO(i , j ) , i  =  l ,nr  ) 

230  continue 
endif 


write e*, 801  l)(p4(i)  ,i=«l,nr) 

8011  formate/'  Prob  of  Disabling  Mechanical  Damage  (P4)' ,5x,3(f 11.4)) 


write (*, 8ul2)epbei ) ,i=l ,nr ) 

8012  formate/'  Prob  of  Fire  Due  to  tlit  on  Ammunition  ( Pb ) '  ,  3x ,  3  (  f  1 1 . 4  ) ) 


writee*,3013)(piiei , 1) ,  i  *  l,nr  ) 

8013  formate/'  Probability  of  Perforating  Sinks  ePll)  Sink  l',3(fll.4)) 
if(nsinks  .eq.  2)  writee*,8014)  e  pll(i,2),  i«l,nr  ) 

8014  format(  '  Sink  2',3efll.4)) 


ife  nr  .gt.  1  )  then 
wri tee*, 8015 ) 

8015  formate/'  Probability  of  Igniting  Fluid  Spilled  on  Set  F', 

*  '  of  Previous  Bursts,'/ 

*  '  Given  a  Perforation  of  tne  Appropriate  Sink  [P12eF)]') 
writee*,801b) 

8016  formate/35x, 'F  -  tl',6x,'F  -  t2',5x,'F  *  tlt2'/) 
do  250  k  *  l.nsintts 


8017 

do  240  i  »  2, nr 

m  -  2**ei-l)  -  1 

if (i.eq.2)writee*,8017)i,k, (pl2ek,i,j ) ,j-l,m) 
format ( 14x, ' Burst ' ,i2, ' ,  Sink ' ,i 2, 5x ,f 7 .4, 2( 5x , ' 

--  ')) 

8018 

if ei.eq.3)write(*,8018)i ,k,(pl2ek,i,j),j«l,m) 
formate  lax, 'Burst' ,i2, ' ,  Sink' ,i2,3(5x,f7.4)) 

240 

continue 

250 

continue 

-  - 

endif 

write(*,8019)  today 

37 


8019  format (//’  This  run  was  made  '  ,aou) 


do  270  j  =  1,6 

do  2ou  i  *  1,9 
pCj  ,i )  -  0. 

260  continue 

270  continue 

write(*,8128)  identl 
write(*,8002)  idenc2 
write(*,8021) 

8021  format(/29x, '0  UTPUT  LIST',/) 
do  290  j  -  1,28 
do  280  k.  -  1,6 
s(j  ,fc)  »  0.0 
280  continue 

290  continue 

c 

if (icor .ne . 1 )then 

c  perform  calculations  for  case  of  DEPENDENT  events 
c 

c  combine  p9's  and  plO's  appropriately, 
c 

do  320  k  =  1,  2 
do  310  i  *  1,  nr 
pi(i,k)  *  0. 
p5(i,k)  -  0. 

if(  k.  .eq.  1  .or.  nsinKs  .eq.  2  )  then 
do  300  j  ■  1,  nis(k) 
m  *  isink.(K,  j  ) 
pl(i  ,k)  *  p  1  ( i  ,1c)  +  p9(i,m) 
p5(i,k.)  =  p3(i,K.)  +  p9(i,m)  *  pl0(i,m) 

300  continue 

endif 

310  continue 
320  continue 


c 

c  ******  calculate  state  probabilities  for  first  burst  ****** 

C  *★★*★****★*  *  it  it -kit  irk  it*  ***  DEPENDENT  *************  ★**★★**★***★** 

c 

call  pdl 

c----------------------------  -- 

if(  nr  .gt.  1  )  then 
c 

c  ******  calculate  state  probabilities  for  second  burst  ****** 
c  it************* *************ijepenDENT****** ****************** 

c 

c  -----  —  -  --  --  --  --  --  --  ----------- 

call  pd2 

c----------  -  ------------------- 

endif 

if (  nr  .gt.  2  )  then 
c 

c  *****  calculate  state  probabilities  for  thira  burst  ***** 
c  *************** ******** DEPENDENT**** ********************* 
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call  pd3 


endif 

write (*, 8022) 

8022  f ormat (// 19x, 'Calculations  Made  Assuming  Dependent  Events'// 

*  42x , ' N  UMBER  OF  BURST  S'// 

*  47*,'l',9x,'2',9*,'3'/> 

write  (*,8023)  (p(i , 1) ,i*l ,nr )  - 

8023  format(8x, ' 1.  Undamaged  ' ,3(fl0.5)) 

write  (*,8024)  (p(i ,2) ,i*l ,nr ) 

8024  format(8x, ' 2.  Puncture  Only  ',3(fl0.5)) 

write  (*,8025)  (p(i ,3) ,i*l,nr) 

8025  f ormat ( 8x , ' 3.  Mecnanical  Only  ’,3(fl0.5)) 

write  (*,8028)  ( pv i , 4 ) , i= i ,nr ) 

8G2o  f ormat ( 8x ,' *♦.  Fire  umy  ',3(fl0.5)) 

i£(  nr  .eq.  1  )write  (*.8027) 

8027  f ormat ( 8x , ' 5.  Puncture,  Mechanical,  No  Fire','  —  ') 

if(  nr  .ne.  1  )write  (*,8028)  (p(i ,5) ,i*2,nr ) 

8U28  f ormat ( 8x, ' 5.  Puncture,  Mechanical,  No  Fire','  —  ', 

*  2(f 10. 5) ) 

if(  nr  .eq.  1  )write  (*,8029) 

8029  £ormat( 8x, ' 6.  Mechanical,  Fire  ','  —  ') 

if(  nr  .ne.  1  )write  (*,8030)  (p(i ,6) ,i»2,nr ) 

8030  format (8x, '6.  Mechanical,  Fire  —  ', 

*  2(f 10. 5) ) 

write  (*,8UJi)  (p ( i , 7 ) , i-1 , nr ) 

8031  £ornac(8x, ' 7.  Mechanical  ’,3(fl0.5)) 

write  (*,8U32)  ( p( i , 8 ) , i= 1 , nr ) 

8032  f ormat ( 8x ,  '  8.  Fire  ' ,3(fl0.5)) 

write  (*,80j3)  (p(i , 9) , i»l  ,nr  ) 

8033  format(8x, '9.  Either  Mechanical  or  Fire  ',3(fl0.5)) 

do  330  i-l,nr 
pbrl ( i  ,  1)  -  p(i  ,  1) 

pbrl (i ,2)  «  p ( i , 2 )  +  p(i,3)  +  p(i,5) 
pbrl (i  , 3)  »  p(i,9)  +  p(i,2.) 
pbrivi  ,4)  =»  p(i  ,8) 

330  continue 

write(*,8034)  (pbrl (i , 1) ,i*l ,nr ) 

8034  format(/9x, 'Neither  M,  nor  F,  nor  K  (1)  ',3(fl0.5)) 

write(*,80o5)  (pbrl (i , 2) ,i-l  ,nr ) 

8035  format(9x,’M  or  F,  but  not  K  (2  +  3  +  5) ' ,3(f 10.5) ) 
write(*,803b)  (pbrl (i , 3) , i- 1 , nr ) 

8036  format(9x,'M  or  F  (including  K)  (2  +  9)  ’,3(fl0.5)) 

write(*  ,8037 )  (pbrl (i  ,  4 ) , i<*  1  ,nr ) 

8037  format (9x,'K  (8)  ' ,3(f 10.5)//) 

c 

endif 

ir(,icor.ne.0)tnen 

c  perform  calculations  for  case  of  independent  events 
c 

do  350  j  =*  1,3 

do  3^*0  i  *  1,9 
p( j ,i  )  -  0. 
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J4U  continue 

350  continue 

c  combine  p9's  and  plO’s  appropriately, 

c 

do  380  k  -  1,  2 
do  370  i  ■  1,  nr 
ql(i,k)  -  1. 
q5(i,k)  -  1. 

if(  k  .eq.  1  .or.  nsinks  .eq.  2  )  then 
do  360  j  »  1,  nis(k) 
m  *  isink(k , 1 ) 

ql(i,k)  -  ql(i.k)  *  (  1.  -  p9(i,m)  ) 
q5(i,k)  -  q5(i,k)  *  (  1.  -  p9(i,m)  *  plG(i.m)  ) 
360  continue 

endif 

pl(i,it)  -  1.  -  ql(i,k) 
p  5  ( i  ,  k )  *  i.  -  q5(i,k) 

37u  continue 
3b0  continue 
c 

c  ******  calculate  state  probabilities  for  first  burst  ****** 
c  ***********************  independent**** ********************* 


call  pil 


if  (  nr  .gt.  1  )  then 
c 

c  ******  calculate  state  probabilities  for  second  burst  ****** 
£  *************************  INDEPENDENT********************1**** 


call  pi  2 


endif 

if  (  nr  .gt.  2  )  then 
c 

c  *****  calculate  state  probabilities  for  tnird  burst  ***** 

c  **********************  independent************************ 


call  pi3 


endif 


write(*,8038) 

80o8  format(// ldx, 'Calculations  Made  Assuming  Independent  tvents'// 

*  42x, 'NUMBER  OF  BURSTS'// 

*  *+7x,’l',9x,'2',9x,'3'/) 
write  (*,8023)  (p(i , 1 ) ,i-l,nr ) 
write  (*,3024)  (p( i  , 2 ) , i»  1 , nr ) 
write  (*,8025)  (p(i , 3) ,i*l,nr ) 
write  (*,8026)  (p(i ,4) ,i-l,nr ) 
write  (*,8039)  (p(i ,5) ,i-l,nr) 

8039  format ( ax ,' 5.  Puncture,  Mecnanical,  No  Fire ' , 3(f 10. 5) ) 


,3(110.5)) 


write  (*,8040)  (p(i , 6 ) ,i*l ,nr ) 

8U40  format ( 8x ,' o.  Mechanical,  Fire 
write  (,*,8031)  (p( i  ,  7 )  ,i-i ,nr , 
write  (*,8032)  ( p( i , 8 ) , i  =  1 , nr ; 
write  (*,8033)  (p(i ,9) ,i*i,nr ) 
do  390  i*l,nr 
pbrl(i,l)  -  p(i , 1 ) 

pbrl (i ,2)  -  p(i ,2)  +  p(i,3)  +  p(i,5) 
pbrl(i ,3)  -  p(i ,9)  +  p(i,2) 
pbrl(i,4)  *  p(i,8) 

390  continue 

write (* , 8034)  (pbrl (i , 1) ,i«l ,nr ) 
write(*,8035)  (pbrl(i ,2),i*l,nr) 
write(*,803b)  (pbrl (i , 3) ,i«l ,nr ) 
write(*,8037)  (pbrl(i ,4) ,i*l,nr ) 
endif 

write (*, 8019 )  today 

if (ipv.eq.O)goto  410 
do  400  i  “1,2 

if(i.eq.l  .and.  icor .eq . 1 )goto  400 
if(i.eq.2  .and.  icor .eq . 0)goto  400 
ifd.eq.l  .or.  nsinics  .eq.2)then 
write(*,8128)  identl 
write(*, 8002)  ident2 
else 

write(*,8041) 

8041  format (///// ) 
endif 

if(i  .eq.  1)  write (*, 8042) 

8042  format ( lox, ' Individual  State  Probabilities  -  Dependent  Events') 

if(i  .eq.  2)  write(*,8043) 

8C-*_>  format ( 15x Individual  State  Probabilities  -  Independent  Events') 
write (* , 8044) 

804*  format (//21x, ' -  Burst  1  - ',7x,' -  Burst  2  - ',7x, 

*  ' - Burst  3 - ' ) 

write (*, 8045) 

8045  i ormat (/ Ox, 'F  G',lUx,'g  “0  g  *  l',8x,'g  “0  g  *  1' 

*  ,  8x , '  g  »  0  g  -  1 '  / ) 
if (  i  .eq.  1  )  then 

write(*,8046)s(l,l),s(l,4),s(l,2),s(i,5),s(l,3),s(l,6) 

8046  format (5x, 'phi  phi ' ,2x, 3( 4x , 2f 9. 5) ) 

if  (  nsinks  .eq.  2  )  then 

write (* , 8047 )s(33,l),s(33,2),s(33,5),s(33,3),s(33,6) 

8047  format (4x, '  phi  t 1 ' , 7x ,f 9. 5, •  --  ’ , 2(4x, 2f 9.5) ) 

write (* , 8048 )s(9,2),s(9,5),s(9,3),s(9,6) 

8048  format (4x, '  phi  t2',7x,2('  --  ' ) , 2( 4x , 2f 9.5) ) 

write (* , 8049 )s (3,3),s(3,b) 

8049  format( 5x, 'pni  t3' ,3x,2(4x,2( '  —  ' ) ) , 4x, 2f 9.5) 


write (* , 8050)s (4l , 2) ,s( 41 , 3) ,s ( 41 ,o) 


8050  format( 5x, 'phi  tlt2,,5x,2('  .  —  '),4x,f9.5,'  — 

*  4x ,  2f  9 . 5 ; 
c 

write(*,805i)s(35,3),s(35,6; 

8051  format(5x,  'pni  tit3'  ,5x,2( '  —  ’),4x,2C’  —  ’),4x,2f9.5) 

c 

write (* , 8052 )s (ll,3),s(ll,6) 

8052  format(5x, ’pni  t2t3' ,5x,2( '  —  '),4x,2('  —  ’),4x,2f9.5) 

c 

write (*,8053 )s(43, 3) 

8053  format( 5x, 'phi  tlt2t3* ,3x,2(2( ’  —  ' ) ,4x) ,f 9.5, '  -  ’) 

endif 

c 

write(*, 8054)s( 17,l),s(17,2),s(17,5),s(17,3),s(l7,6) 

8054  formate 5x,'  t 1  phi ' ,6x,f 9. 5, '  —  ’ ,2(4x,2f 9.5)) 

c 

if  (  nsinxs  .eq.  2  )  tnen 
write(*,8055)s(25,2),s(25,3) ,s(25,6) 

8055  format (bx, ' 1 1  t2',7x,2e'  —  ’ ),4x,f9.5,'  —  '  ,4x,2f9.5) 

c 

write (*, 8056)s(19, 3 ),s( 19,6) 

8056  format (bx, 'tl  t3',7x,2('  —  *),4x,2(’  —  ’),4x,2f9.5) 

c 

write(*,8u57)s^27,3) 

805/  f ormat ( 6x , ' 1 1  1 2t  3 ' ,  5x ,  2(  2( '  —  ' ) ,4x),f 9.5, '  —  ') 

endif 
c 

write<,*,8058)s(5,2),s(5,5),s(5,3),s(5,6) 

8058  formate 5x,’  t2  phi’, 6x,2('  —  ’ ) , 2(4x ,2f 9.5) ) 

c 

if (  nsinxs  .eq.  2  )  then 
writee*,8059)s(37,2),se37,3),s(37,6/ 

8059  f ormat ( bx , ' t 2  tl',7x,2(’  —  ’),4x,f9.5,’  —  ’,4x,2f9.5) 

c 

write (* , 8060 )s( 7 , 3) ,s ( 7 , 6) 

80b0  f ormat ( bx , ' t 2  t3',7x,2(’  —  *),4x,2('  —  ’),4x,2f9.5) 

c 

write (* , 8061 )s ( 39 , 3) 

80ol  f ormat ( bx t2  tlt3' ,5x,2(2( '  —  * ) ,4x) ,f 9. 5, ’  —  ') 

endif 
c 

write(*,b0b2)s(2,3) ,s(2,6) 

8062  format(5x,'  t3  phi',6x,2('  —  '),4x,2('  — 

*  2(  4X ,  2f  9, 5) ) 
c 

if  (  nsinks  .eq.  2  )  then 
write(*,80b3)s(34,3),3(34,6) 

8063  f ormat ( bx , ' t 3  tl',7x,2('  —  ’),4x,2('  —  ' ),4x,2f9.5) 

c 

write(*,8064)s( 10, 3) ,s( 10,6) 

806h  format (ox, 't3  t2',7x,2('  —  '),4x,2('  —  ';,4x,2t9.5 

c 

write ( * , 8065 )s ( 42, 3) 

8065  formate  ox, 't3  t It  2 ' , 5x , 2( 2( '  —  ' ) ,4x)  ,r 9. 5 , '  — 

endif 


42 


wri  ce  ( * ,  80oo  )s  ( 2  i ,  2 ) ,  s  ( 21 , 3 ) ,  s  ( 2 1 , 6 ) 

806b  format(4x, 't lt2  pni’,6x,2('  —  ' ) , 4x ,f 9. 5, '  — 

*  -*x ,  2f  9 .  5 ) 

c 

if  (  nsinks  .eq.  2  )  cnen 
write(*,80b7)s(23,3) 

8U67  format(4x, 't lt2  t3*,7x,2(2('  —  ' ) ,4x) ,f 9. 5, '  —  * 

endif 
c 

write(*,b068)s(  18,3)  ,s(  18,6)  -  — 

8068  format(4x, 'tlt3  phi',6x,2('  —  '),4x,2('  —  '), 

*  4x,2f 9.5) 
c 

if  (  nsinks  .eq.  2  )  then 
write(*,80b9)s(26,3) 

bOo*  format ( 4x, ' tit  3  c2',7x,2(2('  —  '  )  ,  4x) ,f 9. 5, '  —  ' 

endif 
c 

write (*,  8070 )s(  6,  3)  ,s  (.6, 6; 

807'J  format  ( 4x, 't2t  3  pni',6x,2('  —  '  ),4x,2('  —  ’), 

*  4x ,  2f  9. 5 ) 
c 

if  (  nsintcs  .eq.  2  )  then 
write (*,807 l;s(3b,3) 

8071  format(4x, 't2t3  tl',7x,2(2('  —  ' ) , 4x) ,f 9. 5, '  —  ' 

endif 

c 

wri te (*, 8072 )s(  22,3) 

8072  format(2x, 'tit2t3  phi ' ,6x,2(2( '  —  ' ) , 4x) ,f 9. 5, '  — 

c 

write  (.* ,  8u73  )s  ( 65, 1 )  ,s  ( t>5, 2)  ,s  ( 65, 5)  ,s  ( 65, 3 )  ,s  ( 65, 6) 

8073  f ormat ( 8x, '* ' , 9x,f 9. 5, ( '  —  ' ) , 2( 4x, 2f 9.5) ) 

c 

if  (  nsinks  .eq.  2  )  write (*, 8019 )  today 
else 
c 

c  start  here  for  INDEPENDENT  events 

c 

write (* , 8046 )t (l,l),t(l,4),t(l,2),t(l,5),t(l,3),t(l,6) 
c8046  format( 5x, 'phi  phi ' , lx , 3( 4x , 2f 9. 5) ) 
c 

if  (  nsinks  .eq.  2  )  then 

write(*,8074)t(33, l),t( 33,4), t( 33,2), t( 33,5), t( 33,3), t(  33, 6) 

8074  format(4x,'  phi  t 1 ' , 7x , 2f 9. 5, 2( 4x , 2f 9 .5) ) 
c 

wri te (* , 8048 )t ( 9, 2) ,t(9,5),t(9,3),t(9,6) 
c8048  format(4x,'  phi  t2',7x,2('  —  ' ) , 2( 4x , 2f 9.5) ) 

c 

write ( * , 8049 )t ( 3, 3) , t ( 3, 6) 

c8U4y  format ( 5x ,' phi  t 3 ' , 2x, 2( 4x , 2( '  —  ' ) ) ,4x,2f 9.5; 

write(*,8075)t (41 ,2) , t ( -1, 5) , t(-l , 3; , t(4i ,6) 

8075  f ormati 5x, 'pni  tlt2',5x,2^'  —  ’ ) ,2(4x,2f 9.5) ) 


write(*,8051)t(35,3),t(35,6) 

c8051  f ormat ( 5x , 'phi  t  lc3 '  ,  5x ,  2(. '  —  *),4x,2('  —  '), 

c  *  4x ,  if  9 . 5 ) 

c 

write(* ,8052 )t ( li,3),t(il,o) 

c8052  format(  5x, 'phi  t2t3',5x,2('  —  '),4x,2('  —  ' ) ,' 

c  *  4x ,  2f  9 . 5 ) 

c 

write(*,807b)t(43,3),t(43,6) 

8076  format(5x, ’phi  t It2t3 ' , 3x , 2( 2( '  —  ' ) ,4x) ,2f 9.5) 

endif 

c 

write(*,8077)t(17,l),t(17,4),t(17,2),t(17,5),t(17,3),t(17,6) 

8077  format( 5x, '  tl  phi ' , 6x, 2f 9. 5 ,2e4x ,2f 9. 5) ) 
c 

if  (  nsinks  .eq.  2  )  tnen 

write  (*,8078)te49,I),t(49,4),t(49,2),t(49,5),t(49,3j,t(49,6) 

8078  f ormat (5x,'  tl  t i ' , 7x , 2f 9. 5 , 2( 4x , 2f 9. 5) ) 
c 

writee*,8U79)te25,2),te25,5),te25,3),te25,6) 

8079  formate6x, 'tl  t2' ,7x,2e '  —  ' ) , 2e 4x , 2f 9.5) ) 

c 

write(*, 805b)t ( 19,3),t(19,6) 

c805o  format^6x , 't 1  t3',7x,2('  —  '),4x,2('  —  ’),4x,2f9.5) 

c 

write(*,8080)t( 57 ,2),t(57,5),te57,3),te57,6) 

8080  formate6x, 'tl  tlt2',5x,2e'  —  * ) ,2(4x,2£ 9.5) ) 

c 

write e* , 8081 )te  51 , 3) ,t (51,6) 

8081  f ormat ( 6x , ' t 1  tlt3' ,5x,2(2( '  —  ’ ) ,4x) ,2f 9.5) 

c 

write (*, 8082 )t( 27, 3) , t(27,6) 

8082  format( 6x, 't 1  t2t3' ,5x,2(2( '  —  ' ) , 4x) , 2f 9.5) 

c 

write (* , 8083 )t ( 59, 3), t( 59, 6) 

8083  formate 6x, 't 1  t  It2t3 ' , 3x , 2( 2( '  —  ’ ) ,4x) ,2f 9. 5) 

endif 

c 

write (* ,8U58)t( 5, 2) ,t(5,5),te5,3),t(5,6) 
c8058  format(5x,'  t2  pni',6x,2('  —  ' ) ,  2( 4x  ,  2f 9.5) ) 

c 

if  (  nsinks  .eq.  2  )  then 

write (* , 8084)t ( 37,2),te37,5),t(37,3),t(37,6) 

8084  format(6x, 't2  tl',7x,2e'  —  ' ) ,2(4x,2f 9.5)) 

c 

write (*, 8085 )t( 13,2),t(13,5),t(13,3),t(13,6) 

8085  format ( 6x, 't2  t2',7x,2e'  —  ' ) ,2(4x,2f 9.5) ) 

c 

write (*,8060 )t ( 7,3),t(7,6) 

c8060  format ( 6x , 't 2  t3*,7x,2('  —  '),4x,2('  —  '),4x,2f9.5, 

c 

write(*, 8086)t ( 45 , 2) , t ( 45, 5) , t ( 45, 3) , t (45, 6) 

8086  f ormat ( bx ,' t2  tlt2',5x,2('  —  '  )  ,  2( 4x  ,  2f 9.5) ) 


c 


write ( * , 8087 )t(39,3),t(39,6) 


8087  formate bx, 't2  t  It 8 ' , 5x , 2( 2( ' 


—  ' ) , 4x) , 2f  9. 5 ) 

write  t  * , 8088 )t ( 15,3),t(15,o) 

8uc8  formatC  bx, 't2  t2t3',5x,2e'  —  '  ),4x,2('  —  ' ) , **x , 2r  9. 5) 

write (* , 8089 )t ( 47 , 3) , t (*7 , 6) 

8089  formatC 6x, 't2  t It 2t3 ' , 3x , 2( '  --  ’),4x,2(’  —  ’) 

*  ,4x, 2f  9. 5) 

endif 

c 

writee*. 8062)t ( 2, 3) , t  ( 2, 6)  - 

c8062  format(5x,'  t3  phi',6x,2e'  —  '),4x,2('  '), 

c  *  2(4x,2f 9.5)) 

c 

if  (  nsinks  .eq.  2  )  then 
writee* , 8063 )t ( 34 , 3) , t (  34 , 6) 

c80b3  f ormat e 6x , ' t 3  tl',7x,2('  —  '),4x,2('  '),4x,2f9.5) 

c 

write  e*,8064)t(10,3),t(10,6) 

cbOo*  f ormat e ox , ' t 3  t2',7x,2e'  —  ’),4x,2e’  —  ’),4x,2f9.5) 

c 

write (* , 8090)t e 4 , 3 ) , t ( 4, 6) 

8090  format e Ox , ' c 3  t3',7x,2('  —  '),4x,2('  ’),4x,2f9.5) 

c 

writee*. 8091)t(42,3),t(42,6) 

8091  formatehx, 't3  tlt2' ,5x,2(2( ’  —  ' ) ,4x) ,2f 9.5) 

c 

wri te ( * , 6092 ) t ( 3b , 3 ) , t ( 36 , 6 ) 

8092  format ( bx, 't3  tlt3',5x,2(2e'  —  ' ) ,4x) ,2f 9.5) 

c 

writee*. 8093 )t ( 12,3),t( 12,6) 

8U93  formatebx, 't3  t 2t 3 ' , 5x , 2( 2( '  —  ' ) , 4x ) , 2f 9 . 5 ) 

c 

writee* , 8094)t  e44, 3),te44,6) 

8094  f ormat e 6x , ' t 3  t It2t3' ,3x,2(2( '  —  ' ) ,4x) ,2f 9.5) 

endif 

c 

writee*,8095)t(21,2),te21,5),te21,3),te21,b) 

8095  f ormat ( 4x, ' t It 2  pni'>6x,2('  —  * ) , 2( 4x,2f 9.5) ) 

c 

if  (  nsinks  .eq.  2  )  then 

writee*,8096)te53,2),t(53,5),te53,3),te53,6) 

809b  f ormat ( 4x , ' t It 2  tl',7x,2e'  —  ' ) , 2e 4x, 2f 9. 5) ) 

c 

write(*,8097)te29,2),te29,5),te29,3),t(29,b) 

8097  formate4x, 't lt2  t2',7x,2('  —  ' ) ,2(4x,2f 9.5)) 

c 

write (*, 8098 )t ( 23 , 3) , t  e  23, 6) 

8098  f ormat (4x, 'tit 2  t3',7x,2e2e'  —  ' ) ,4x) ,2f 9.5) 

c 

write  e*,8099)tebl,2),tebl,5),tebl,3),tebl,b) 

8099  f ormat ( 4x, ' t It 2  tlt2,,5x,2('  —  ' ) ,2e^x,2f 9.5) ) 

c 

writee* ,81 00)te  55,3),t(55,6) 

8100  formate ^x, 't lt2  1 1 1 3 ' , 5x , 2( 2( ' 


' ) ,  4x ) ,  2f  9 . 5 ) 


~  *  *  *  "V*  ~  *  -  «  -W"  ---  r  r  .'.  v.  . 


write  (*,  61  31,6) 

8101  format (‘♦x,  'tit 2  tzt3'  ,5x,2(2( ' 

write!'  ,  8I02)t  !6_>,  3; , t  (63, 6 ) 

8102  formatOx, 'tlt2  Clt2t3'  ,3x,2(2( ' 


'),4x),2r9.5) 


' ) ,  4x ) ,  2f  9 . 5 ) 


write!*, 8103) 

8103  f orma t ( / / / 5 lx , ' Table  Continued  On  Next  Page'/) 
write!* ,8019)  today 

write!*, 8128)  ident 1 
write!*, 8002)  ident2 
write!*, 8043) 
write!*, 8104) 

8104  f ormat ( 35x ,'( continued )' ) 
write!*  ,8044 ) 

write!*, 6045) 


wri te !*, 8066 ;t ! 18,3), t! 18,6) 

c806d  format !4x, 'tlt3  phi',6x,2!' 
c  *  4x , 2f  9 . 5 ) 


'  ) » 4x ,  2! ' 


it  (  nsinxs  .eq.  2  )  tnen 

write ( * , 8105)t ( 50, 3) , t ( 50, 6) 

81 J5  f ormat ! 4x , ' t It 3  tl',7x,2('  —  ’),4X,2!’  —  ’),4x,2f9.5) 


write (* , d 106 )t ( 26 , 3) , t( 26 , 6 J 

8106  format ! 4x , ' t It 3  t2',7x,2(2!’ 

write !*, 8107 )t! 20, 3), tv  20, 6) 

81 07  format(4x, 'tlt3  t3',7x,2!2(' 

write!* , 8108 )t ! 56, 3) , t ! 58, 6) 

8108  format! 4x , ' t It 3  t lt2' ,5x,2( 2( ' 

write ! * , 8109)t ( 52 , 3 ) , t (52, o) 
6109  format ( 4x , ' 1 1 1 3  t it 3 ' , 5x , 2(2( ’ 


—  ’ ) ,  4x ) ,  2f  9 . 5 ) 


' ) ,  4x ) ,  2f  9 . 5 ) 


'  ) ,  4x ) ,  2f  9 . 5 ) 


’)  ,4x;,2f 9.5) 


write!*,dli0)t(28,3),t!28,6) 

8110  format! 4x, 't It 3  t2t3' ,5x,2(2( 1  —  ' ) ,4x) ,2f 9.5) 


write(*,8111)t(60,3),t(60,6) 

8111  format (4x, 'tit 3  tlt2t3' ,3x,2(2( ' 
end  if 


' ) , 4x ) , 2f  9 . 5 ) 


write (*, 8070)t ( 6, 3) , t ( 6.6) 

c8070  format(4x, 't2t3  phi',6x,2(' 
c  *  4x ,  2f  9. 5 ) 


,),4x,2(' 


if  !  nsinxs  .eq.  2  )  tnen 

write!*,6112)tl38,3),t(38,6) 

8112  t ormat ( 4X , ' t 2t 3  tl,,7x,2!2l'  —  ' ) ,4x) ,2f 9. 5) 

write(*,8113)t( 14, 3) , t ( 14, 6) 
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* ) ,  4x ) ,  2f  9 . 5 ) 


8113  f oraat ( 4x , ' t 2t 3  t2',7x,2(2(' 

write  (*  ,  81 1-t  )t  (  8, 3  ) ,t  (8, 6) 

8114  format ( 4x, 't2t 3  t3',7x,2(2(' 

write (*, 81 15)t (46,3) ,t (4b, 6) 

8115  format (4x, 't2t3  t lt2 ' , 5x , 2( 2( ' 

write(*,8116)t(40,3),t(40,6) 

8110  format ( 4x t2t3  tlt3' ,5x,2(2( • 

write(*,8117)t(16,3),t(16,6) 

8117  format(4x, 't2t3  t 2t 3 1 , 5x , 2( 2( ' 

write (*, 81 18)t ( 48, 3) , t ( 48 , 6) 

8118  format(4x, 't2t3  t It 2t 3 ' , 3x , 2( 2( ' 
endif 


’ ),4x),2f9.5) 


—  ' ) ,  4x ) ,  2f  9 . 5 ) 


—  ' ) ,  4x ) ,  2f  9 . 5 ) 


—  ' ) ,  4x ) ,  2f  9 . 5 ) 

—  ' ) ,  4x ) ,  2f  9 . 5 ) 


c 

wri te (* , 81 19;t(22,3),t(22,b) 

8119  format(2x, 'tlt2t3  phi ' ,6x,2(2( * 
c 

if  (  nsinXs  .  eq.  2  )  then 

wri te (* , 8 12u)t ( 54 , 3) , t ( 54 , 6) 

8120  format ( 2x t lt2t 3  ti',7x,2(2(' 
c 

write  (*,  812 1  )t(30,o),t(.30,6) 

8121  format ( 2x, 't It2t3  t2',7x,2(2(' 
c 

write(*,8122)t(24,3),t(24,6) 

8122  f oraat ( 2x , ' t lt2t 3  t3',7x,2(2(' 
c 

write (* , 8123)t ( o2, 31 , t ( 62, o) 

8123  format ( 2x , ' t It 2t  3  t lt2 ’ , 5x , 2( 2(  ’ 
c 

wri te( * ,8124)t(56,3),t(56,b) 

8124  format ( 2x , ' t it 2t 3  t It 3 ' , 5x , 2( 2( ' 
c 

write(*,8125)t(J2,3),t(J2,b) 

8125  format ( 2x , ' t it 2t 3  t2t3',5x,2l2(.’ 


—  ' ) ,  4x ) ,  2f  9 . 5 ) 


—  ' ) ,  4x ) ,  2f  9 . 5 ) 


—  ' ) ,  4x ) ,  2f  9 . 5 ) 


—  '  ) ,  4x ) ,  2f  9 . 5 ) 

—  ' )  ,4x)  ,2f 9.5) 


) ,  4x ) ,  2f  9 . 5 ) 


—  ' ) ,  4x ) ,  2f  9. 5 ) 


write (* , 812b)t ( 64, 3) , t (b4, 6) 

8126  format ( 2x, 't lt2t 3  t It 2t3 * , 3x , 2( 2( '  ~  ' ) , 4x) ,2f 9.5 ) 

endif 


c 

write (*, 8127 )t(b5,l),t(65,4),t(65,2),t(65,5),t(65,3),t(65,6) 
8127  format ( 8x , , 5x , 3(4x,2f 9.5) ) 
c 

write(*,8019)  today 
endif 

400  continue 

4  i  G  continue 

stop 
end 
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subroutine  pd 1 
c 

c  ******  calculate  state  probabilities  for  first  burst  ****** 

£  ★★*******************L)£p£^££^‘x  EVENTS  ******* *************** 

c 

implicit  real*8  (a-h,  o-z) 

common  pl(3,2),  p3(2,2:3,3),  p4(3),  p5(3,2),  p(6,9),  p6(3), 

*  P9( 3, 10) ,  p 10(3, 10) ,  pll(3,2) ,  pl2(2,2: 3,3) , 

*  q 1 ( 3, 2 ) ,  q3( 2, 2 : 3, 3 ) ,  q4(3),  q5(3,2),  q6(3), 

*  8(65,6),  t(65,6),  nis(2),  isink(2,10)  - 

c  phi  phi  0  Table  5  Formula  1 

c 

s(l,l)  -  1.  -  pl(  1, 1)  -  p 1 ( 1 , 2 )  -  p4{ 1 )  -  p6( 1 ) 
c 

P(l,l)  =  p(l,l)  +  s(l,  1) 

c  tl  pm  0  Table  5  Formula  2 

c 

s ( 1 7 ,  i )  =  pl(l,l)  -  p5(  1 , 1 ) 
c 

P(l,2)  =  p ( 1,2)  +  s ( 1 7 , 1 ) 

c  -  --  --  --  --  --  --  --  --  --  --  --  --  --  --  - 

c  pni  tl  0  Table  5  Formula  3 

c 

s(33,  i)  =  p  1  ( 1 , 2)  -  p5( 1 , 2) 
c 

P( 1 , 2)  »  p( 1,2)  +  s ( 33, 1 ) 

c  -  --  --  --  --  --  --  --  --  --  --  --  --  --  --  - 

c  star  0  Table  5  Formula  4 

c 

s(65,l)  -  p 5 ( 1 ,  1 )  +  p 3 ( 1,2)  ■+•  pb(  1 ) 
c 

P(  1,4)  -  p(  1,4)  -r  s  (  65,  1) 
p C 1,8)  -  p ( 1,8)  +  s ( 65, 1 ) 
p( 1, 9)  -  p(  1 ,9)  +  s( 65, 1 ) 

c  phi  pm  1  Table  5  Formula  5 

c 

s (  1 , 4  )  =  p4( 1 ) 
c 

P (  1,3)  »  p( 1 , 3)  +  s ( 1 ,4) 

P C  1,7)  -  p<.  1,7)  +  s ( 1 , 4) 

P( 1,9)  -  p ( 1,9)  +  s ( 1 , 4 ) 

c  -  --  --  --  --  --  --  --  --  --  --  --  --  --  --  - 
c  tl  phi  1  Table  5  Formula  6 

c 

s(  17,4;  »  0. 
c 

P( 1 , 5)  =  p ( 1,53  +  s ( 17,“) 

Pk  1 , 7 )  =  p  v  1,7)  +  s ( 17 , l ) 

P (  1,9)  =  p ( i , 9 )  +  s( 17,4) 

c  pni  tl  1  Table  5  Formula  7 

c 
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subroutine  pd2 
c 

c  ******  calculate  state  probabilities  for  second  burst  ****** 
c  ** ******** ****** ******* DEPENDENT  EVEiVTS  **************  ******* 

c 

implicit  real*8  (a-'n,  o-z) 

common  pl(3,2),  p3(2,2:3,3),  p4(3),  p5(3,2),  p(6,9),  p6(3), 

*  p9( 3, 10) ,  pl0(3,10),  pll(3,2) ,  pl2( 2,2: 3,3), 

*  q 1 ( 3 , 2 ) ,  q3(2, 2:3,3),  q4(3),  q5(3,2),  q6(3), 

*  s(65,6),  t(65,6),  nis(2),  isinlt(2,10) 

c------------------------------ _ 

c  phi  phi  0  Table  5  Formula  1 

c 

s( 1, 2)  -  8(1,1)  *  (  1.  -  p 1 ( 2 , 1 )  -  p 1 ( 2 , 2 )  -  P4( 2)  -  p6( 2)  ) 
c 

P( 2, 1 )  -  p( 2, 1 )  +  s ( 1 , 2) 


c  t2  pni  0  Table  5  Formula  2 

c 

s ( 5, 2)  =  s(l,l)  *  (  pl(2, 1)  -  p5( 2, 1 )  ) 
c 

p(2,2)  =  P(2,2)  +  s ( 5, 2) 

c  -  --  --  --  --  --  --  --  --  --  --  --  --  --  --  - 
c  phi  t2  0  Table  5  Formula  3 

c 

s( 9, 2)  =  s(l,l)  *  (  pl(2,2)  -  p5( 2, 2)  ) 
c 

p(2,2)  -  p( 2, 2)  +  s ( 9, 2) 

c  tl  phi  0  Table  5  Formula  1 

c 

sU7,2)  »  s ( 17 , 1 )  *  (  1.  -  pi ( 2,  1 )  -  pi 1 2, 2) 

*  -  p 3 ( 1,2,1)  -  P4(2)  -  P6( 2)  ) 
c 

p(2,2)  -  p( 2, 2 )  +  s( 17,2) 

c  tlt2  pni  0  Table  5  Formula  2 

c 

s ( 21 , 2 1  -  s ( 1  7  , 1 )  *  (  p 1 ( 2, 1 )  -  p5(2,  1 )  ) 
c 

p( 2,2)  -  p( 2, 2)  +  s( 21 ,2) 

c  cl  t2  0  Table  5  Formula  3 

c 

8(25,2)  -  s(  17, 1)  *  (  pl( 2, 2)  -  p5( 2, 2)  ) 
c 

P( 2 , 2 )  -  p( 2, 2)  +  s ( 25, 2) 

c  pni  tl  0  Table  5  Formula  1 

c 

s(33, 2)  -  s ( 33, 1 )  *  (  1.  -  p  1  ( 2 , 1 )  -  pl(2,2) 

*  -  p3( 2, 2, 1 )  -  p4(2)  -  P6( 2)  ) 
c 

p( 2, 2)  -  p( 2, 2)  +  s(33,2) 
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phi  tlc2  0 


Table  5  Formula  3 


s(41,2)  -  s(3i, 1)  *  (  pl(2,2)  -  p5(2,2)  ) 


P( 2, 2)  -  p( 2, 2)  +  s(4i,2) 


p(2,5)  *  p(2, 5)  +  3(17,5) 
p(2,7)  -  p(  2, 7)  -h  s  ( 17 , 5 ) 
p(2,9)  -  p( 2, 9)  +  s(  17,5) 


c 

c 

c 

c 


c 

c 

c 

c 


c 
c 
c 

I 

c 


I  c 

c 

c 

c 

c 

c 

c 

I  c 

T 


i 


c 

c 

c 


c 


clc2  pni  1  Table  5  Formula  6 

s(21,5)  -  s( 17,4)  *  (  p 1 ( 2, 1 )  -  P5(2, 1)  ) 

P( 2, 5)  -  p( 2, 5)  +  8(21.5) 

P(2, 7)  -  p(2,7)  +  s(21,5) 

P( 2, 9)  -  p( 2, 9)  +  s(21,5) 


cl  c2  1  Table  5  Formula  7 

s(25,5)  *  s(17,4)  *  (  pl(2,2)  -  p5(2,2)  ) 

P( 2, 5)  =  p( 2, 5)  +  s(25,5) 
p( 2, 7)  =  p( 2, 7)  +  s ( 25, 5) 
p( 2, 9)  -  p( 2, 9)  +  s ( 25, 5) 


phi  cl  1  Table  5  Formula  5 

s(33,5)  =■  s ( 33, 1 )  *  P4( 2)  +  s(33,4)  *  (  1.  -  pl(2,l) 
-  pl(2,2)  -  p3( 2,2, 1 )  -  p6( 2)  ) 

P( 2, 5)  =  p( 2, 5)  +  s ( 33, 5) 

P( 2, 7)  -  p( 2, 7 )  +  s ( 33, 5) 

P( 2, 9)  -  p(2, 9)  +  s ( 33, 5) 


c2  cl  1  Table  5  Formula  6 

s(  37, 5)  =*  s  (  33, 4)  *  (  pl(2,l)  -  p5(2,l)  ) 

p(2,5)  -  p( 2, 5)  +  s(37,5) 

P( 2, 7)  -  p( 2, 7)  +  s(37,5) 

P( 2, 9)  -  p( 2, 9)  +  s ( 37 , 5) 


phi  Clc2  1  Table  5  Formula  7 

s ( 41 , 5)  -  s ( 33 , 4 )  *  (  pl(2,2)  -  p5(2,2)  ) 

p( 2, 5)  -  p(2, 5)  +  s ( 41 , 5 ) 
p( 2, 7)  -  p( 2, 7 )  +  s(41,5) 
p( 2, 9)  -  p( 2, 9)  +  s ( 41 , 5) 


scar  1  Table  5  Formula  6 

s(65,5)  *  s(65,l)  *  p4(2)  +  s(65,4) 

+  s(l,4)  *  (  p5( 2 , 1 )  +  p5( 2,2)  +  p6( 2)  ) 
+  s(17,4)  *  (  p5( 2, 1 )  +  p5( 2, 2) 

+  p3(  1 , 2 , 1 )  +  p6(2)  ) 

+  s ( 33, 4 )  *  (  p5( 2, 1 )  +  p5( 2, 2 ) 

+  p3( 2, 2, 1 )  +  P6(2)  ) 


p( 2, 6)  *  p( 2, 6)  +  s( 65,5) 


2,7;  =  p  (.  2 , 7 )  +  s(6o,5) 
2,8;  -  pC2,a;  +  3(65,5; 
2,9)  =  p( 2, 9)  +  s ( 65 , 5 ) 


ptoc  =  p( 2, 1 )  +  p( 2 , 2)  +  p( 2, 3)  +  p( 2 , 4)  +  p(2,5)  +  p(2,6) 
if(  aos(  i.  -  ptot  )  .gt.  .000001  )  write  (*,100)  ptot 
100  format ('  sum  of  DEPENDENT  probabilities  for  2  bursts  is 
return 


subroutine  pd3 


c 

c  *****  calculate  state  probabilities  for  tnird  Durst  ***** 

C  ********************J}£PENJj£N’£  EVENTS********************* 

C 

implicit  real*8  (a-h,  o-z) 

common  pl(3,2)„  p3(2,2:3,3),  p4(3),  p5(3,2),  p(b,9),  po(3), 

*  p9( 3, 10) ,  P10( 3, 10) ,  pll(3,2),  pl2( 2, 2: 3,3) , 

*  ql(3,2),  q3(2,2:3,3),  q4(3),  q5(3,2),  q6(3), 

*  s( 65,6) ,  t(65,6),  nis(2),  isinlc(  2, 10) 


c  pni  phi  0  Table  5  Formula  1 

c 

s( 1 , 3)  -  s(l,2)  *  (  1.  -  pl(3, 1)  -  p 1( 3, 2)  -  p4( 3)  -  p6(3)  ) 
c 

p(  3, 1 )  <■  p(  3, 1 )  +  s(l,3) 

c  ------------------------------- 

c  t3  phi  0  Table  5  Formula  2 

c 

s(2,3)  -  s( 1,2)  *  (  p  1  ( 3, 1 )  -  p5( 3, 1 )  ) 
c 

P( 3 , 2 )  -  p( 3, 2 )  +  s ( 2, 3) 

c---------  --------------------- 

c  phi  t3  0  TaDle  5  Formula  3 

c 

s(3,3)  *  s(l,2)  *  (  p 1 ( 3, 2 )  -  p5(3, 2)  ) 
c 

P( 3, 2)  -  p( 3, 2)  +  s ( 3, 3) 

c------------------------------ 

c  t2  pm  U  TaDle  5  Formula  1 

c 

s(5,3)  =  s ( 5, 2)  *  (  1.  -  pl( 3, 1)  -  pi ( 3, 2) 

*  -  p3( 1,3,2)  -  p4( 3)  -  P6( 3)  ) 
c 

P( 3, 2)  -  p( 3, 2)  +  s( 5, 3) 

c------------------------------ 

c  t2t3  phi  0  Table  5  Formula  2 

c 

s(o, 3)  »  8(5,2)  *  (  p 1 ( 3, 1 )  -  p5( 3, 1 )  ) 
c 

P( 3, 2 )  -  p( 3, 2)  +  s ( 6, 3 ) 

c------------------------------ 

c  t2  t3  0  Table  5  Formula  3 

c 

s(7,3)  -  s ( 5, 2)  *  (  p 1 ( 3, 2)  -  p5(3,2)  ) 
c 

P( 3, 2)  -  p( 3, 2)  +  s ( 7, 3) 

c------------------------------ 

c  phi  t2  0  TaDle  5  Formula  1 

c 

s(9,3j  -  s(9,2)  *  (  1.  -  p 1 ( 3 ,  1 )  -  p 1 ( 3, 2 ) 

*  -  p3( 2, 3, 2)  -  P4( 3)  -  p6( 3)  ) 
c 

P( 3, 2)  *  p( 3, 2)  +  s ( 9, 3) 

c  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  ^  —  —  —  —  —  —  —  —  — 
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V 


t3  t2  0  Tabie  5  Formula  2 

sU0,3)  =■  s(9,2)  *  (  p  1  ( 3, 1 )  -  p5(3, 1)  ) 

P( 3, 2)  =  p( 3, 2)  +  s (10,3) 

phi  t2t3  0  Taole  5  Formula  3 

s(ll,3)  -  s(9, 2)  *  (  pl(3,2)  -  p5(3,2)  ) 

P( 3, 2)  -  p(3,2)  +  s( 11,3) 

tl  phi  0  Table  5  Formula  1 

s(17,3)  -  s ( 17 , 2)  *  (  1.  -  p 1 ( 3, 1 )  -  p 1 ( 3, 2) 

*  -  p 3 ( 1,3,1)  -  p4( 3)  -  p6( 3)  ) 

P( 3, 2)  *  ?( 3, 2)  +  s ( 17 , 3) 

tit 3  phi  0  Table  5  Formula  2 

s(la,3)  =  s ( 17 ,2)  *  (  pi ( 3, 1 )  -  p5(3, 1)  ) 

p( 3, 2)  =  p( 3, 2)  +  s( 18, 3) 

ci  t3  (J  Taole  5  Formula  3 

s(19,3)  -  s( 17,2)  *  (  pi ( 3, 2)  -  p5(3,2)  ) 

P  (  3 , 2 )  =  p( 3, 2)  +  s( 19, 3) 

tit 2  phi  0  Tabie  5  Formula  i 

s(21,3)  -  s(21,2;  *  (  1.  -  pl( 3, 1 )  -  pl( 3, 2) 

"  -  p3( 1 , 3, 3)  -  p4(3)  -  p6( 3)  ) 

P(3,2)  =  p( 3, 2)  +  s(21,3) 

tlt2t3  phi  0  Table  5  Formula  2 

s ( 22 , 3)  -  s ( 21 , 2)  *  (  p 1 ( 3, 1 )  -  p5(3,l)  ) 

P(3, 2)  -  p ( 3, 2 )  +  s ( 22, 3) 


tit 2  t3  0  Table  5  Formula  3 

s ( 23, 3)  *  s ( 21 , 2 )  *  (  pi C 3, 2)  -  p5(3,2)  ) 
P( 3, 2)  -  p( 3, 2)  +  s ( 23, 3 ) 
tl  t2  0  Table  5  Formula  1 


s( 25, 3)  =  s ( 25, 2)  *  (  1.  -  p 1 ( 3, 1 )  -  pl(3,2)  -  p3(l,3,l) 

-  p3( 2,3,2)  -  p4(3)  -  p6( 3)  ) 


t  It  3  1 2  0 


P(3,2)  +  s ( 25 , 3) 
Table  5  Formula  2 


c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


S<26,3)  *  5(25,2)  *  (  pl(3,l)  -  p5(3,l)  ) 

P( 3, 2)  -  p( 3, 2)  +  s( 26, 3) 
tl  t2t3  0  Table  5  Formula  3 

s(27,3)  -  s(25, 2)  *  (  pl(3, 2)  -  p5(3,2)  ) 

P ( 3,2)  -  p( 3, 2)  +  s ( 27, 3) 

0  Table  5  Formula  1 

s(33,2)  *  (  1.  -  p 1 ( 3, 1 )  -  pl(3, 2) 

-  p3( 2, 3, 1 )  -  P4(3)  -  pb( 3)  ) 

p(  3 , 2 )  =■  p(  3, 2)  +  s  ( 33, 3 ) 
t3  tl  0  Table  5  Formula  2 

s ( 34 , 3)  -  s ( 33, 2)  *  (  pi ( 3, 1 )  -  P5(3,l)  ) 

P( 3 , 2 )  -  p( 3, 2)  +  s ( 34, 3) 
phi  tit 3  0  Table  5  Formula  3 

s(35,3)  =  s ( 33,2)  *  (  p 1 ( 3, 2)  -  P5(3,2)  ) 

P( 3 , 2)  -  p( J, 2)  +  s(35, 3) 
t2  tl  U  Table  5  Formula  1 

s(37,3)  -  5(37,2)  *  (  1.  -  p 1 ( 3, 1 )  -  pi(3,2)  -  p3( 1,3,2) 

-  p3( 2, 3, 1 )  -  P4( 3)  -  pb( 3)  ) 

p ( 3,2)  -  p( 3, 2)  +  5(37,3) 
t2t3  tl  0  Table  5  Formula  2 

6(  38,3)  -  5  (  37,2)  *  (  p  1  ( 3 ,  1)  -  P5(3,l)  ) 

P( 3, 2)  -  p( 3, 2)  +  s(3b,3) 
t2  tit 3  0  Table  5  Formula  3 

6(39,3)  -  5(37,2)  *  (  pl(3,2)  -  p5(3,2)  ) 

P( 3, 2)  -  p( 3,2)  +  s(39,3) 
pni  tit2  0  Table  5  Formula  1 


phi  1 1 

s( 33, 3)  - 


s(4i,3)  -  s ( 4 i , 2 )  *  (  1.  -  pl(3, 1)  -  p  1  ( 3 , 2 ) 


-  p3(  2 , 3 , 3)  -  P4(3)  -  P6( 3 )  ) 


c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


c 


c 

c 

c 


c 


c 

c 

c 

c 


p(  3,2)  =*  p(  3, 2)  +  s(4l,3) 


t3  tit  2  U  Taoie  5  Formula  2 

s(42,3)  -  s ( 41 , 2)  *  (  pl(3, 1)  -  p5(3,l)  ) 
P( 3,2)  -  P( 3, 2)  +  s(42,3) 


phi  tlc2c3  0  Table  5  Formula  3 

s(43,3)  -  s ( 41 , 2 )  *  (  pi(3,2)  -  p5(3,2)  ) 
p(3,2)  -  p( 3, 2)  +  s(43,3) 


star  0  Table  5  Formula  4 

s ( 65 , 3)  =  5(bo,2)  *  q4(3) 

+  s ( 1 , 2)  *  (  p5( 3, 1 )  +  p5( 3, 2)  +  p6l 3 )  ) 
+  s(  17,2)  *  (  p5( 3, 1 )  +  p 5 ( 3,2) 

-i-  p3(  1,3,1)  +  pb(  3)  ) 

'  +  s( 5,2)  *  (  p5( 3,1)  +  P5( 3,2) 

1  +  P3(  1,3,2)  +  p6( 3)  ) 

'  +  s ( 21 , 2)  *  (  P5( 3,1)  +  p5( 3, 2) 

+  p3( 1,3,3)  +  p6( 3)  ) 

+  s(33,2)  *  (  p5( 3,1)  +  p5(3,2) 

+  p3( 2, 3,1)  +  P6( 3)  ) 

+  s(9,2)  *  (  P5( 3, 1 )  +  p 5 ( 3,2) 

+  P3(2, 3,2)  +  p6( 3)  ) 

+  s(41,2)  *  (  P5( 3, 1 )  +  P5(3,2) 

+  p3( 2, 3, 3)  +  p6( 3)  ) 

+  s(25,2)  *  (  p5( 3, 1 )  +  p5( 3,2) 

-r  p  3  (  1,3,1)  +  p3(  2, 3, 2)  +  p6(3)  ) 

+  s( 37 , 2)  *  (  P5( 3, 1)  +  p5( 3,2) 

^  p3( 1,3,2)  +  p3( 2, 3,1)  +  p6( 3)  ) 

P( 3 , 4 )  -  p( 3, 4)  +  s( 65, 3) 
p(3,8)  =■  p(3,8)  +  s(65,3) 

P( 3,9)  -  p( 3, 9)  +  s ( 65, 3) 


phi  phi  1  Table  5  Formula  5 

8(1,6)  *  8(1,2)  *  p4( 3)  +  s( 1,5)  *  (  1.  -  pl( 3, 1 ) 
-  pl( 3,2)  -  P6( 3 )  ) 

p(3,3)  *  p( 3, 3)  +  s( 1,6) 
p( 3, 7)  -  p ( 3 , 7 )  +  s ( 1 , 6) 

P( 3, 9 )  -  p( 3, 9)  +  s( 1,6) 


t3  pni  1  Table  5  Formula  6 

s( 2, o)  -  s( 1, 5)  *  (  pit  3, 1 )  -  p5(3, 1 )  ) 


P( 3, 5)  *  p( 3, 5)  +  s( 2, 6) 


p( 3,7)  =  p( 3 , 7 )  + 
P(3,9)  *  p( 3, 9)  -r 


1 

C 

pni  t3  1  Taole  5 

Formula  7 

c 

s(2,6)  =*  3(1,5)  *  (  p  1  ( 3, 2) 

-  p5( J, 2)  ) 

> 

»> 

o 

£ 

c 

P(  3, 5)  -  p( 3, 5) 

+  s(3,6) 

P( 3, 7)  -  p(3,7) 

+  s(3,6) 

p( 3, 9)  -  p( 3, 9) 

+  s ( 3, 6) 

t2  pni  L 


Table  5  Formula  5 


s(5,6)  -  s ( 5, 2)  *  P4(3)  +  s(5, 5)  *  (  1.  -  pi(3,l) 
*  -  p 1(3,2)  -  p3( 1 , 3, 2)  -  p6( 3)  ) 

P(  3, 5)  -  p( 3, 5)  +  s ( 5, 6) 
p(  3, 7)  -  p(  3, 7)  4-  s(5,6) 

P( 3, 9)  =  p( 3, 9)  +  s ( 5, 6) 

t2t3  pni  1  Table  5  Formula  6 

s  ( 6 , 6 )  >  s ( 5, 5 )  *  (  p  L( 3, 1 )  -  p5(3,i)  ) 

p( 3, 5 )  -  p( 3, 5)  +  s(6, 6) 
p( 3, 7)  -  p( 3, 7)  +  s(6,6) 

P( 3, 9)  -  p( 3, 9)  +  s ( 6, 6 ) 


t2  t3  1  Table  5  Formula  7 

s(7,6)  =  s(5,5)  *  (  P 1  ( 3 , 2 )  -  P5( 3,2)  ) 

p( 3,5)  *  p(3,5)  +  s ( 7, 6 ) 
p(3,7)  -  p( 3, 7 )  +  s(7,6) 
p( 3, 9)  -  p(3,9)  +  s ( 7, 6) 


pni  t2  1 


Table  5  Formula  5 


s(9, 6)  =  s(9,2)  *  P4(3)  +  s ( 9, 5)  *  (  1.  -  pl(3,l) 
*  -  pl( 3,2)  -  p3(2,3,2)  -  P6( 3)  ) 


c 

c 


tl  t2  1 


Table  5  Formula  5 


9(25,6)  -  s(25,2)  *  p4(3)  +  s(25,5)  *  (  1.  -  pl(3,l) 

*  -  p 1 ( 3 , 2)  -  p3( 1,3,1)  -  p3( 2 , 3, 2)  -  p6(3)  ) 
c 

P( 3, 5)  -  p( 3, 5)  +  8(25,6) 
p(3,7)  -  p(3,7)  +  8(25,6) 

P( 3, 9)  -  p( 3, 9)  +  8(25,6) 

c  clt3  t2  1  Table  5  Formula  6  ~ 

c 

s(26,6)  -  s(25,5)  *  (  pl(3,l)  -  p5(3,l)  ) 
c 

p( 3, 5 )  -  p( 3, 5)  +  8(26,6) 
p(3,7)  =  p ( 3 , 7 )  +  s(2b,6) 
p( 3, 9)  -  p( 3, 9)  +  s(2b,b) 

c  -  --  --  --  --  --  --  --  --  --  --  --  --  --  --  - 
c  £ i  c2c3  1  Table  5  Formula  7 

c 

s(27,b)  -  a(25,5)  *  (  p 1 ( 3 , 2 )  -  p5(3,2)  ) 
c 

P(3, 5)  =  p( 3, 5)  +  s(27,6) 

P(3,7)  «  p( 3,7)  +  8(27,6) 

P( 3, 9)  =  p( 3, 9)  +  s(27,6) 

c------------------------------ 

c  phi  cl  1  Table  5  Formula  5 

c 

s(33,6)  =  s(33,2)  *  p4(3)  +  s(33,5)  *  (  1.  -  pl( 3, 1 ) 

*  -  pl( 3,2)  -  p3(2,3,l)  -  p6( 3)  ) 
c 

p( 3, 5)  =  p( 3, 5)  +  s ( 33 , 6) 
p( 3, 7)  *  p( 3, 7 )  +  s(33,6) 

P( 3, 9)  »  p( 3, 9)  +  s(33,6) 

c  c3  cl  1  Table  5  Formula  6 

c 

9(34,6)  =■  8(33,5)  *  (  p 1 ( 3 , 1 )  -  p5(3,l)  ) 
c 

P( 3, 5)  -  p( 3, 5)  +  s(34,6) 

P(3,7)  -  p(3,7)  +  s(34,6) 

P(3,9)  -  p( 3 , 9)  +  s(34,6) 

c  phi  clc3  1  Table  5  Formula  7 

c 

s(35,6)  -  s ( 33, 5)  *  (  p 1 ( 3, 2)  -  p5(J,2)  ) 
c 

P( 3, 5)  -  p(3, 5)  +  8(35,6) 

P(3,7)  -  p( 3, 7)  +  s(35,6) 

P( 3,9)  -  p( 3, 9)  +  s(35,6) 

c-_---------------------------- 

c  £2  cl  1  Table  5  Formula  5 

c 

s (  37 , 6 J  =  s ( 3/ ,  2)  *  p4(3)  -r  s(37,5)  *  (  1.  -  pl(3,l) 

*  -  pl( 3, 2)  -  p3( 1,3,2)  -  P3( 2, 3, 1 )  -  p6(3)  ) 


fll] 


c2t  3  cl  1 


P(3»  5)  =*  p(  3, 5)  +  s  (  3 
p( 3 , 7 )  -  p( 3 , 7 )  +  s  (  3 
p(3,9)  =  p ( 3 , 9 )  +  s ( 3 

Taole  5  Formula  6 


s  (  38 , 6 )  -  s(37,5)  *  (  p 1 ( 3, 1 )  -  P5( 3, 1 )  ) 

P( 3, 5)  -  p( 3, 5)  +  s ( 38 , 6) 
p(3,7)  -  p( 3, 7)  +  s(38,6) 

P( 3, 9)  -  p( 3, 9)  +  s ( 38, 6 ) 


c  lt3  1 


Table  5  Formula  7 


s ( 39 , 6)  =  s( 37,5)  *  (  pl(3,2)  -  p5(3,2)  ) 

P( 3, 5)  =  p(3, 5)  +  s ( 39 , 6) 
p( 3 , 7 )  -  p( 3, 7)  +  s ( 39 , 6 ) 

P( 3, 9)  =  p( 3, 9)  +  s ( 39 , 6) 


phi  Clt2  1 


Table  5  Formula  5 


s(41,6)  -  s ( 41 , 2 )  *  p4( 3)  +  s ( 41 , 5 )  *  (  1.  -  pi  ( 3,  1 ) 
'  -  p  1( 3, 2)  -  p3(2,3,3)  -  P6(3)  ) 

p( 3, 5)  -  p( 3, 5)  f  s ( 41 , 6) 
p(3 , 7 )  »  p(3,7)  +  s(41 ,6) 
p ( 3,9)  -  p( 3, 9)  +  s ( 41 ,6) 


cj  t  It  2  1 


Table  5  Formula  6 


s  ( >42 , 6 )  =  s(4l,5)  *  (  p  1  ( 3, 1 )  -  p5(  3, 1 )  ) 

P( 3 , 5 )  -  p( 3 , 5 )  +  s(42,6) 
Pt3,7)  -  p( 3, 7)  +  s(42,6) 
pU,9)  -  p (  3,9)  +  s  ( 42 ,  b ) 


phi  tlc2t3  1 


Table  5  Formula  7 


s(43,o)  -  3(41,5)  *  (  p 1 ( 3 , 2 )  -  P5( 3, 2)  ) 

p( 3, 5 )  »  p( 3,5)  +  s(4j,6) 
p ( 3 , 7 )  -  p( 3 , 7 )  +  s ( 43 , 6 ) 
p( 3, 9)  -  p( 3, 9)  +  s(43,6) 


star  1 


Table  5  Formula  8 


s(65,b)  *  s(65,2)  *  p4(3)  +  s(65,5) 

+  s(l, 5)  *  (  p5( 3, 1 )  +  p5( 3 , 2 )  +  P6(3)  ) 
t  s(  17,5)  *  (  p5(  3, 1 )  -4-  p5( 3, 2) 

+  p3( 1 ,3, 1)  +  p 6 i 3 )  ) 

*  s (,  5 , 3 )  *  (  p5(  3,  1 )  t  ?5(3,2) 

+  p3( 1,3,2)  +  po( 3 )  ) 

+  s  (  21 , 5)  *  (  p5(  3,  1 )  +  p5(.3,2) 
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subroutine  pi  1 


c  ******  calculate  state  probabilities  for  first  burst  ****** 

c  **************x******xnjdEPENDENT  events******************** 

c 

implicit  real*8  (a-n,  o-z ) 

common  pl(3,2),  p3(2, 2:3,3),  p4(3),  p5(3,2),  p(6,9),  p6(3), 

*  P9(3,10),  plO( 3,10),  pll(3,2) ,  pl2(2,2: 3,3), 

*  q 1 ( 3, 2 ) ,  q3( 2, 2: 3,3),  q4(3),  q5(3,2),  q6(3), 

*  s(b5,6),  t(65,6),  nis(2),  isinlc(2,10) 


q  lq  2  -  ql(l,l)  *  q 1 ( 1 , 2)  *  q  6 (  1 ) 

plq2  -  (  pl( 1, 1)  -  p 5 ( 1,1)  )  *  qi( 1 ,2)  *  q6(  1 ) 

P2q 1  »  (  p  1  ( 1,2)  -  p5(  1,2)  )  *  ql(  1 , 1)  *  q6(l) 

plp2  -  (  pl(l,l)  -  P5( 1 , 1 )  )  *  (  pl( 1 , 2)  -  P5( 1,2)  )  *  q 6 ( 1 ) 


c  pni  phi  0  Table  6  Formula  1 

c 

t(i,  1)  =  qlq2  *  q4( 1) 

P(  l,  1 )  -  p(i,l)  +  tci.l) 


tl  phi  0  Table  6  Formula  2 

t ( 17 , 1 )  =  plq2  *  q4(l) 

p(l,2)  *  p(l,2)  +  t( 17, 1) 


pni  tl  0  Table  6  Formula  3 

t ( 33 , 1 )  =  p2q 1  *  q 4 ( 1 ) 

P( 1.2)  =  p ( 1,2)  +  t ( 33, 1 ) 


tl  tl  0  Taole  6  Formula  4 

t ( 49 , 1  )  =  p ip2  *  q 4 (  1 ) 

PC 1,2)  -  p(  1 , 2 )  +  t ( 49 , 1 ) 


*  0  Table  6  Formula  5 

t(65, 1 )  -  q 4 C 1 )  *  (  1.  -  q5( 1 , 1 )  *  q5(l,2)  *  q6(l)  ) 

P(  1,4)  *  p(l,4)  +  t(65, 1) 
p( 1,8)  -  p( 1 ,8)  +  t( 65, 1 ) 
p(  1,9)  *  p ( 1,9)  +  t(65, 1) 


pni  pni  1  Table  6  Formula  6 

t(i,u)  =  qlq2  *  p4( 1 ) 

P C 1,3)  -  p(  1 , 3)  +  t ( 1 , 4 ) 
pi  1,7)  =  p(  1 ,  7)  +  t(  1 ,4) 
p(  1 ,9)  =  pi  1 , 9)  +  t( 1,4) 


tl  phi  I  Table  6  Formula  7 

t ( 17 , 4)  =  piq2  *  p “ (  l  ) 

p(  1, 5)  =*  p(  1,5)  +  c(  17 ,4) 
P(l,7)  -  p(l,7)  +  t C 1 7 ,4) 
p(l,9)  -  p( 1, 9)  +  t( 17,4) 


phi  tl  1  Table  6  Formula  8 

t ( 33, 4 )  =  p2q 1  *  p4( 1 ) 

P(l,5)  -  p( 1 , 5)  +  t(33,4) 

P(  1 ,7)  -  p(  1 , 7 )  +  t( 33,4) 

P ( 1,9)  -  p( 1 , 9)  +  t( 33,4) 


tl  tl  1  Table  6  Formula  9 

t ( 49, 4)  =  p lp2  *  p4( 1 ) 

P ( 1,5)  «  p( 1 , 5)  +  t ( 49 , 4) 

P( 1 , 7 )  =  p(  1 , 7 )  +  t(49,4) 

P( 1, 9)  =  p( 1 , 9)  +  t ( 49 , 4 ) 


*  1  Table  6  Formula  10 

t C 65 , 4 )  *  p4( 1 )  *  (  1.  -  q5( 1 , 1 )  *  q5(l,2)  *  q6(l)  ) 

p(  1,6)  «  p( 1,6)  +  t(65,4) 

P ( 1,7)  -  p ( 1,7)  +  t(65,4) 

P( 1,8)  =  P ( 1 , 8 )  +  t ( 65 , 4 ) 

p( 1,9)  “  p( 1,9)  +  t( 65,4) 


ptot  -  p ( 1 , 1 )  +  p( 1 ,2)  +  p ( 1,3)  +  p( 1 ,4)  +  p( 1,5)  +  p(  1 , 6 ) 
it  (  abs(  1.  -  ptot  )  .gt.  .000001  )write (* , 100)ptot 
100  format( '  sum  of  INDEPENDENT  probabilities  for  1  burst  is  ' 
return 
end 


subroutine  pi2 
c 

c  ******  calculate  state  prooabilities  for  second  burst  ****** 

c  ********************** independent  events******************** 

c 

implicit  real*8  (a-h,  o-z ) 

common  pl(3,2),  p3(2,2:3,3),  p4(3),  p5(3,2),  p(6,9),  pb(3), 

*  p9( 3, 10) ,  p 10(3, 10) ,  pll(3,2) ,  p 12( 2, 2: 3,3) , 

*  ql(3,2),  q3(2, 2: 3,3),  q4(3),  q5(3,2),  q6(3), 

*  s(65,6),  t(65,6),  nis(2),  isink( 2, 10) 

qlq2  -  q 1 ( 2 , 1 )  *  ql(2,2)  *  q6(2) 

plq2  *  (  pl(2, 1)  -  p5( 2, 1 )  )  *  q  1  (  2 , 2 )  *  q6(2) 

P2q 1  -  (  pl(2,2)  -  p5(2,2)  )  *  ql(2,l)  *  q6(2) 

plP2  -  (  pl(2, 1)  -  p5( 2, 1 )  )  *  (  pl(2,2)  -  p5( 2, 2)  )  *  q6(2) 

—  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  — 

c  phi  phi  0  TaDle  6  Formula  1 

t ( 1 ,2)  =  t(l,i)  *  q lq2  *  q4(2) 
c 

P(2, 1)  -  p(2, 1)  +  t( 1,2) 


c  t2  pni  0  Table  6  Formula  2 

c 

t  C  5, 2)  =  t  C 1, 1)  *  Plq2  *  q4( 2) 
c 

P(2,2)  -  p(2,2)  +  t ( 5, 2) 

c  phi  t2  0  Table  6  Formula  3 

c 

t(9,2)  =  t(l,l)  *  p2ql  *  q4( 2) 
c 

P( 2, 2)  =  p(2, 2)  +  t ( 9, 2) 

c  t2  t2  0  Table  o  Formula  4 

c 

t ( 13 , 2)  =  t(l,l)  *  plp2  *  q4(2) 
c 

p( 2,2)  »  p( 2, 2 )  +  t ( 13, 2) 

c------------------------ 

c  tl  pni  0  Table  6  Formula  1 

c 

t( 17,2)  -  t ( 17  , 1 )  *  q lq2  *  q3(l,2,l)  *  q4(2) 
c 

p(2,2)  -  p(2, 2)  +  t ( 17 , 2) 

c------------------------ 

c  tit 2  phi  0  Table  b  Formula  2 

c 

t(21 ,2)  -  t ( 1 7  , 1 )  *  p lq  2  *  q3( 1 , 2, 1 )  *  q4(2) 
c 

P( 2 , 2 )  =  p( 2 , 2 )  +  t ( 21 , 2) 


c  tl  t2  u  Table  b  Formula  3 

c 


t ( 25 , 2)  -  t(17,i)  *  p2q 1  *  q3(l,2,l)  *  q4(2) 


c 


PU.2)  -  p(2,2)  +  t(25,2) 


c  -  —  —  -  —  —  -  —  —  —  _  —  — 

c  tlt2  t2  0  Table  b  Formula  4 

c 

t ( 29 , 2)  -  t ( 17 , 1 )  *  plp2  *  q3( 1 , 2, 1)  *  q4(2) 
c 

P( 2, 2)  -  p(2,2)  +  t(29, 2) 

c  phi  tl  0  Table  6  Formula  1 

c 

t ( 33 , 2)  -  t(33,l)  *  q lq2  *  q3(2,2,l)  *  q4(2) 
c 

P(2,2)  -  p( 2, 2)  +  t ( 33, 2) 

C”----1---  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  ______ 

c  t2  tl  0  Table  6  Formula  2 

c 

c(37,2)  »  t ( 33 , 1 )  *  Plq2  *  q3(2,2,l)  *  q4(2) 
c 

p(2,2)  *  p(2,2)  +  t ( 37 , 2) 

c  —  —  —  —  —  _  —  _  —  —  —  —  —  —  —  —  _  —  _  — 

c  pni  t lt2  0  Table  6  Formula  3 

c 

t(4l,2)  -  t(33,l)  *  P2q 1  *  q3(2, 2,1)  *  q4(2) 
c 

P( 2, 2 J  =  p(2,2)  +  t (41 ,2) 

c  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  _________ 

c  t2  tlt2  0  Table  6  Formula  4 

c 

t(45,2l  =■  t <33, 1)  *  plp2  *  q3(2,2, 1)  *  q4(2) 
c 

P(2, 2)  =  p(  2, 2 )  -r  t  (  45,  2) 

c  ti  cl  0  Table  6  Formula  1 

c 

t(49,2)  -  t ( 49 , 1 )  *  q lq 2  *  q3(l,2,l)  *  q3(2,2,l)  *  q4(2) 
c 

P(  2, 2)  -  p(2,2)  +  t  ( 49  *  2 ) 

c  tlc2  tl  0  Taole  6  Formula  2 

c 

t(53,2)  -  t ( 49 , 1 )  *  plq2  *  q3'l,2,l)  *  q3(2,2,l)  *  q4(2) 
c 

P(2,2)  -  p( 2, 2)  +  t ( 53, 2) 

c  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  _________ 

c  tl  tlt2  0  Table  6  Formula  3 

c 

t( 57,2)  -  t ( 49 , 1 )  *  P2q 1  *  q3(l,2,l)  *  q3(2,2,l)  *  q4(2) 
c 

p(2,2)  -  p( 2, 2)  +  t( 57,2) 
c  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  — 

c  1 1 1 2  t lt2  0  Table  6  Formula  4 

c 

t(0l,2)  -  t(49,l)  *  p lp2  *  q3(  1 , 2, 1 )  *  q3(2,2,l)  *  q4(2) 
c 


P( 2 , 2)  -  p(2,2)  +  t(61,2) 


*  0  Table  6  Formula  5 

t(o5,2)  *  q4(2)  *  (  C( 65, 1 ) 

+  t(l,l)  *  (  1.  -  q5( 2,1)  *  q5(2,2)  *  q6(2)  ) 

+  t ( 17 , 1 )  *  (  1.  -  q5( 2, 1 )  *  q5( 2, 2)  *  q3(l,2,l)  *  q6(2)  ) 

+  t(33, 1)  *  (  1.  -  q5(2, 1)  *  q5(2,2)  *  q3(2,2,i)  *  q6(2)  ) 

+  t(49,i)  *  (  1.  -  q5(2, 1)  *  q5(2, 2) 

*  q3( 1 , 2 , 1 )  *  q3(2,2, 1)  *  q6(2)  )  ) 

P( 2, 4 )  -  p( 2, 4 )  +  t( 65, 2) 
p(z,8)  -  p( 2 , 8 )  +  t ( 65 , 2) 
p(2,9)  -  p( 2, 9 )  +  t(65,2) 


ptii  pni  1  Table  6  Formula  6 

t(l,5)  =  q lq 2  *  (  t  ( 1 ,  1 )  *  p4(2)  +  t(l,4)  ) 

P( 2, 3)  =  p( 2, 3)  +  t( 1,5) 
p(2,7)  =  p( 2, 7)  +  t(  1 ,5) 
p( 2,9)  -  p( 2, 9)  +  t ( 1 , 5) 


t2  pni  1  Table  6  Formula  7 

t(5,5)  -  Plq2  *  (  t(i,l)  *  P4( 2)  +  t ( 1 , 4)  ) 

p(2,5)  =■  p( 2, 5)  +  t( 5, 5) 

P(2, 7)  -  p( 2, 7)  +  t( 5,5) 

P( 2, 9)  -  p( 2,9)  +  t ( 5, 5) 


phi  t2  1  Table  6  Formula  8 

C(9,5)  =  P2q 1  *  (  t ( 1 , 1 )  *  p4( 2)  +  t(l,4)  ) 

P( 2, 5)  =  p( 2, 5)  +  t(9,5) 

P( 2, 7 )  -  p( 2 , 7)  +  t(9,5) 

P( 2, 9)  -  p( 2, 9)  +  t(9, 5) 


t2  t2  1  Table  6  Formula  9 

t( 13, 5)  -  p 1P2  *  (  t(i,l)  *  P4( 2)  +  t ( 1 , 4 )  ) 

P( 2, 5)  *  p( 2, 5 )  +  t ( 13, 5 ) 

P( 2 , 7 )  -  p( 2, 7 )  +  t( 13,5) 

P( 2,9)  -  p( 2, 9)  +  t( 13, 5) 


tl  phi  1  Taole  6  Formula  6 

C ( 17 , 5)  -  q lq2  *  q3(l,2,l)  *  (  t(17,l)  *  p4(2)  +  t(17,4)  ) 

P( 2, 5)  -  p( 2, 5)  +  c C 17, 5) 

P( 2, 7)  -  p( 2 , 7 )  +  t(  17,5) 

P( 2, 9)  -  p( 2, 9 )  +  t( 17,5) 


c 

c 

c 


c 

c 

c 

c 


c 

c 

c 

c 


c 

c 

c 

c 


c 

c 

c 


I  c 

c 
c 

c 

I 


c 

c 

c 


c 


tlc2  phi  1  Table  6  Formula  7 

C(21,5)  =>  plq2  *  q3(  1 ,2, 1 )  *  (  c(17,i)  *  p4(2)  +  t(17,4)  ) 

P( 2, 5)  -  p(2, 5)  +  t(21,5) 
p(2,7)  -  p(2,7)  -c  t( 21 , 5) 

P( 2,9)  -  p(2,9)  -t-  C ( 21 , 5) 


cl  t2  1  Table  6  Formula  8 

t(25,5)  -  P2q  1  *  q3(  1 , 2, 1 )  *  (  t(17,l)  *  p4(2)  +  t(17,4J7 

P( 2, 5)  -  p(2, 5)  +  t ( 25 , 5) 

P( 2, 7)  «  p(2, 7)  +  t(25,5) 
p( 2, 9)  -  p(2, 9)  +  c ( 25, 5) 


clt2  c2  1  Table  6  Formula  9 

c ( 29 , 5 )  =  p  1P2  *  q3(  1 , 2 , 1)  *  (  c(17,l)  *  P4(2)  -r  t(17,4)  ) 

P( 2, 5 )  =  p(2, 5)  +  t ( 29 , 5) 

P( 2 , 7 )  -  P( 2 , 7 )  +  t(29,5) 

P( 2, 9)  =»  p( 2, 9)  +  t(29,5) 


phi  tl  1  Table  6  Formula  6 

t(33, 5)  -  qlq2  *  q3(2,2, 1)  *  (  c(33,l)  *  P4(2)  +  C(33,4)  ) 

p(2,5)  -  p(2, 5)  +  c ( 33 , 5 ) 

P(  2, 7)  -  p(  2, 7)  ■*  c(33,5) 

P( 2, 9)  -  p( 2, 9 )  +  t ( 33, 5) 


C2  t 1  1  Table  6  Formula  7 

t( 37 , 5)  =  Plq2  *  q3(2,2,l)  *  (  t(33,l)  *  p4(2)  +  t(33,4)  ) 

P(2,5)  -  p( 2, 5)  +  t ( 37, 5) 

P(2,7)  -  p( 2, 7)  +  t ( 37 , 5) 

P( 2, 9)  -  p(2,9)  +  t(37,5) 


phi  Clt2  1  Table  6  Formula  8 

t ( 41 , 5)  -  p2q 1  *  q3( 2, 2, 1 )  *  (  t(33,l)  *  p4(2)  +  t(33,4)  ) 

P( 2, 5 )  -  p(2,5)  +  t(41,5) 

P(2,7)  -  p(2,7)  +  t(41 ,5) 

P(2,9)  -  p(2,9)  +  c(41,5) 


t2  tlt2  1  Table  6  Formula  9 

t ( 45, 5)  -  p 1 p 2  *  q3( 2, 2, 1 )  *  (  t(33,l)  *  p4(2)  +  t(33,4)  ) 

P( 2, 5)  *  p( 2, 5 )  f  c(45, 5) 

P( 2, 7)  -  p( 2, 7)  -r  t(45,5) 
p(2,9)  -  p( 2, 9 )  +  t (45, 5) 


I 
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c 

c 

c 


c 


c 

c 

c 


c 


c 

c 

c 


c 


e 

c 

c 


c 


c 

c 


c 


c 


c 


cl  cl  1  Table  6  Formula  6 

t ( 49, 5;  -  qlql  *  q3(l,2,l)  *  q3(2,2,l) 

*  (  t ( 49 , 1 )  *  P4( 2)  +  t ( 49 , 4 )  ) 

p(2,5)  -  p(2, 5)  +  t (49, 5) 
p(2,7)  -  p(2,7)  +  t(49,5) 

P( 2, 9)  -  p(2,9)  +  t(49,5) 


tit 2  tl  1  Table  6  Formula  7  - 

t( 53, 5)  -  plq2  *  q3( 1 , 2 , 1 )  *  q3(2,2,l) 

*  *  (  t(49, 1)  *  P4(2)  +  t(49,4)  ) 

p( 2, 5)  -  p( 2, 5)  +  t ( 53 , 5) 
p( 2 , 7 )  *  p( 2 , 7)  +  t( 53, 5) 
p(  2, 9)  =»  p<  2,9 )  +  t  (  53, 5 ) 


tl  tlti  1  Table  6  Formula  8 

t ( 57 , 5)  *  p2q 1  *  q J(  1 , 2, 1 )  *  q3(2,2,l) 

*  *  (  t(49, 1)  *  p4(2)  +  t ( 49 , 4 )  ) 

P( 2, 5 )  -  p( 2, 5)  +  t ( 57 , 5) 
p(2,7)  -  p( 2 , 7 )  +  t( 57,5) 

P(2,9)  -  p( 2, 9)  +  t ( 57 , 5) 


tlt2  tlt2  1  Table  6  Formula  9 

c ( 61 , 5)  -  p lp2  *  q 3 ( 1 , 2 , 1)  *  q3( 2, 2, 1 ) 

*  *  (  t ( 49 , 1 )  *  P4( 2)  +  t ( 49 , 4)  ) 

p( 2, 5)  -  p( 2, 5)  +  t(61,5) 
p( 2 , 7)  -  p( 2 , 7)  +  t(61 ,5) 
p( 2 , 9 )  -  p( 2, 9)  +  t ( bl , 5) 


★ 


Table  6  Formula  1U 


t ( 65, 5) 

* 

*  *  ( 

* 

*  *  ( 

* 

★  *  ( 

* 


*  ( 


■  t(65, 1)  *  p4(2)  +  t(65,4) 

+  (  t(l,l)  *  P4( 2)  +  t ( 1,4)  ) 
1.  -  q 5( 2, 1 )  *  q5(2, 2)  *  q6(2)  ) 

+  (  t( 17, 1)  *  p4(2)  +  t( 17,4) 
1.  -  q 5( 2, 1 )  *  q5( 2, 2 )  *  q3( 1,2,1)  * 
+  (  t ( 33, 1 )  *  P4( 2)  +  t(33,4) 
1.  -  q 5( 2, 1 )  *  q 5( 2, 2 )  *  q3(2,2,l)  * 
+  (  t (49 , 1 )  *  p4( 2 )  +  t(49,4) 
1.  -  q5( 2, 1 )  *  q5( 2, 2)  *  q3(l,2,l)  * 


) 

q6( 2)  ) 

) 

q6(2)  ) 

) 

q3( 2,2,1)  *  q6(2)  ) 


p( 2, 6 )  -  p( 2, 6 )  +  t( 65, 5) 
p(  2 , 7 )  -  p(. 2 , 7 )  +  t(65,5) 
p( 2 , 8)  -  p( 2, 8)  +  t ( 65, 5) 
pi  2 , 9)  -  p( 2, 9)  +  t(65,5) 


ptot  -  p ( 2 , 1 )  +  p(2,2)  +  p(2,3)  +  p( 2 , 4 )  +  p(2,5)  +  p(2,6) 


I 
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subroutine  pi3 


*****  calculate  state  probabilities  for  third  burst  ***** 

***************** *** independent  events******************* 

implicit  real*8  (a-h,  o-z) 

common  pl(3,2),  p3(2, 2:3,3),  p4(3),  p5(3,2),  p(6,9),  p6(3), 

*  p9(3, 10) ,  p!0( 3, 10) ,  pll( 3, 2) ,  pl2( 2,2: 3,3) , 

*  ql ( 3, 2) ,  q3(2,2:3,3),  q4(3),  q5(3,2),  q6(3), 

*  s(65,6),  t(65,6),  nis(2),  isink(2,10) 


q lq2  -  ql( 3, 1)  *  ql(3,2)  *  q6(3) 

plq2  =•  (  pl( 3, 1 )  -  P5( 3, 1 )  )  *  ql(  3, 2)  *  q6(3) 

P2q 1  -  (  pl( 3,2)  -  p5(3,2)  )  *  ql(3,l)  *  q6(3) 

plp2  -  (  pi ( 3, 1 )  -  P5( 3, 1 )  )  *  (  pl(3, 2)  -  p5(3, 2)  )  *  q6(3) 

q5q6  *  q5(3,l)  *  q5(3,2)  *  q6(3) 

phi  phi  0  Table  6  Formula  1 

t(  1 , 3)  *  t ( 1 , 2)  *  q lq2  *  q4(3) 

P(3, 1)  -  p(3,l)  +  t( 1,3) 

t3  phi  0  Table  6  Formula  2 

t(2,3)  -  t( 1,2)  *  plq2  *  q4( 3) 

p(3,2)  -  p( 3, 2)  +  t( 2, 3) 

pni  t3  0  Table  6  Formula  3 

t ( 3 , 3)  =  t ( 1 ,2)  *  p2q 1  *  q4( 3) 

P( 3, 2)  -  p( 3, 2)  +  t( 3,3) 

t3  t3  0  Table  6  Formula  4 

t(4,3)  «  t( 1,2)  *  pip2  *  q4( 3 ) 

P( 3, 2)  -  p(3,2)  +  t ( 4, 3) 

t2  phi  0  Table  6  Formula  1 

t(5,3)  *  t( 5, 2)  *  qlq2  *  q4(3)  *  q3(l,3,2) 

p(3,2)  -  p(3, 2)  +  t ( 5, 3 ) 

t2t3  phi  0  Table  6  Formula  2 

C ( 6 , 3 )  -  t( 5, 2)  *  plq2  *  q4(3)  *  q3( 1 , 3, 2) 

P( 3, 2)  -  p( 3, 2)  +  t ( 6, 3) 


1 2  1 3  0 


Table  6 


Formula  3 


t3  c2t3  0 


Table  6  Formula  4 


t  ( 12 , 3 )  =  c ( 9, 2)  *  p lp2  *  q4(3)  *  q3(2,3,2) 


P( 3, 2;  »  p( 3,2)  +  c(12,3) 


1 2  t2  0 


Table  6  Formula  1 


t ( 13, 3)  -  r ( 13,2)  *  q lq 2  *  q4(3)  *  q3(l,3,2)  *  q3(2,3,2) 
p( 3, 2)  -  p( 3, 2)  +  c C 13,3) 


t2t  3  c  2  0 


Table  6  Formula  2 


t(l4,3)  -  t(  13,2)  *  plq 2  *  q4( 3)  *  q3(l,3,2)  *  q3( 2, 3,2) 
P( 3, 2)  -  p(3, 2)  +  t( 14,3) 


1 2  1 2t  3  0 


Table  6  Formula  3 


t  ( 15 , 3)  -  t  ( 13, 2)  *  P2q 1  *  q4(3)  *  q3(l,3,2)  *  q3(2,3,2) 
p( 3,2)  =  p( 3, 2 )  +  t( 15, 3) 


c2t  3  t2t3  0 


Table  6  Formula  4 


tC 16,3)  -  c ( 13,2)  *  p lp2  *  q4( 3 )  *  q3(l,3,2)  *  q3(2,3,2) 


tl  pni  0 


Table  6 


Formula  1 


t( 17,3)  -  C ( 17 , 2)  *  q lq 2  *  q4(3)  *  q3(l,3,l) 
p(3,2)  -  p(3,2)  +  t ( 17 , 3) 
tic 3  pni  0  Table  b  Formula  2 

t ( 13, 3)  -  t ( 17, 2)  *  plq2  *  q4(3)  *  q3(l,3,l) 

P( 3, 2)  -  p( 3, 2)  +  t ( 18 , 3 ) 
tl  t3  0  Table  6  Formula  3 

t(ly,3)  =  t(  17,2)  *  p2q 1  *  q4( 3)  *  q3( 1 , 3 , 1 ) 

P( 3, 2)  -  p(3, 2)  +  t ( 19 , 3) 
tit 3  t3  0  Table  6  Formula  4 

t ( 20, 3 )  =  t ( 17 , 2)  *  p lp2  *  q4( 3)  *  q3(l,3,l) 

P( 3,2)  -  p( 3, 2)  +  t ( 20, 3) 
tlt2  pni  0  Table  b  Formula  1 

t ( 21 , 3)  =  t ( 21 , 2)  *  q lq2  *  q4(3)  *  q3(l,3,3) 

P( 3 , 2 )  -  p( 3, 2)  +  t( 21 , 3) 
tlt2t3  phi  0  Table  6  Formula  2 

t  C  22 , 3 )  =  t ( 21 , 2)  *  p lq  2  *  q4(3)  *  q3(l,3,3) 

PC 3, 2)  «  p( 3,2)  +  t ( 22 , 3 ) 
tlt2  t3  0  Table  6  Formula  3 

t ( 23 , 3)  -  t (21 , 2)  *  P2q 1  *  q4(3)  *  q3(l,3,3) 

P( 3, 2)  -  p( 3, 2)  +  t(23, 3) 
tlt2t3  t3  0  Table  6  Formula  4 

c ( 24 , 3 )  -  t ( 21 , 2)  *  p lp2  *  q4( 3)  *  q3( 1,3,3) 

P( 3, 2)  -  p( 3, 2 )  +  t ( 24 , 3 ) 
tl  t2  U  Table  6  Formula  1 

t ( 25 , 3 )  =  t ( 25 , 2 )  *  q lq 2  *  q4(3)  *  q3( 1 , 3, 1 )  *  q3(2,3,2) 


t lt3  t2  0 


Table  6 


Formula  2 


c 

t(26,3)  -  t(25, 2)  *  piq2  *  q4(3)  *  q3(l,3,l)  *  q3(2,3,2) 
c 

p( 3,2)  -  p( 3, 2)  +  t(26,3) 

c  tl  t2t3  0  Table  6  Formula  3 

c 

t(27,3)  -  t(25,2)  *  P2q 1  *  q4(3)  *  q3(l,3,l)  *  q3(2,3,2> 
c 

P(3,2)  «  p(3,2)  +  t(27,3) 

c  tlt3  t2t3  0  Table  6  Formula  4 

c 

t(28,3)  =  t(25,2)  *  p lp2  *  q4(3)  *  q3(l,3,i)  *  q3(2,3,2) 
c 

p( 3, 2)  -  p( 3, 2)  +  t ( 28 , 3 ) 

c  Clt2  c2  0  Taole  6  Formula  1 

c 

t(29,3)  =  t(29,2)  *  qlq2  *  q4(3)  *  q3(l,3,3)  *  q3(2,3,2) 
c 

p(3,2)  =  p(3,2)  +  t( 29, 3) 

c  clt2t3  t2  0  Table  6  Formula  2 

c 

t(30,3)  =*  t(29,2)  *  plq2  *  q4( 3)  *  q3(l,3,3)  *  q3(2,3,2) 
c 

P( 3,2)  -  p(3,2)  +  t( 30, 3) 

c  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  — 

c  tlt2  t2t3  0  Table  6  Formula  3 

c 

t(  31 , 3)  **  t  ( 29 , 2)  *  P2q  1  *  q4(3)  *  q3(l,3,3)  *  q3(2,3,2) 
c 

p( 3, 2)  =  p(3,2)  +  t(31, 3) 

c  tlt2t3  t2t3  0  Table  b  Formula  4 

c 

t(32,3)  -  C ( 29 , 2)  *  plp2  *  q4(3)  *  q3(l,3,3)  *  q3(2,3,2) 
c 

p( 3, 2 )  -  p( 3, 2)  +  t ( 32, 3) 


c  phi  tl  U  Table  6  Formula  1 

c 

t ( 3o , 3)  -  t ( 33, 2)  *  q lq2  *  q4(3)  *  q3(2,3,l) 
c 

p( 3,2)  -  p( 3, 2)  +  t ( 33, 3) 
c  ________________________ 

c  t3  cl  u  Table  b  Formula  2 

c 

t ( 3a , 3 )  -  1 1 33, 2 )  *  p Iq 2  *  q4(3)  *  q3(2,3,l) 
c 


p( 3 , 2 )  -  p( 3, 2)  *  t ( 34 , 3 ) 


pni  ticJ  0  Table  6  Formula  3 

t ( 35, 3)  =  t(3o,2)  *  p2q 1  *  q4(3)  *  q3(2,3,l) 
p( 3,2)  -  p( 3, 2)  +  t(35, 3) 
t3  tlt3  0  Table  6  Formula  4 

t ( 36 , 3 )  *  t ( 33,2)  *  p 1P2  *  q4(3)  *  q3(2,3,l) 

P( 3, 2 )  -  p ( 3 , 2)  +  t(36,3) 
t2  cl  0  Table  6  Formula  1 

C ( 37 , 3 )  -  C ( 37 , 2)  *  q lq2  *  q4(3)  *  q3(l,3,2)  *  q3(2, 3, 1 ) 
p(  3, 2)  =*  p(  3, 2;  +  c  ( 37, 3) 
t2t3  tl  J  Table  6  Formula  2 

t(  38, 3)  =  c ( 37 , 2)  *  plq2  *  q<*(3>  *  q3(l,3,2)  *  q3(2,3,l) 
p( 3, 2)  -  p(3,2)  +  t( 38, 3) 
t2  clc3  0  Table  6  Formula  3 

c(39,3)  =  C ( 37 ,2)  *  P2q 1  *  q4(3)  *  q3(l,3,2)  *  q3(2,3,l) 
p( 3, 2)  -  p( 3 , 2)  +  t ( 39 , 3) 
c2c3  clc3  U  Table  6  Formula  4 

t(40,3)  =  t (37,23  *  p lp2  *  q4( 3)  *  q3(l,3,2)  *  q3(2,3,l) 
p( 3, 2)  -  p  (  3 , 2 )  +  t ( 40, 3 ) 
phi  tlc2  0  Table  6  Formula  1 

t(4i , 3)  >  t( 41 ,2)  *  qlq2  *  q4(3)  *  q3(2,3,3) 

P(3,2)  -  p( 3,2)  +  t ( 41 , 3) 
t3  tlt2  0  Table  6  Formula  2 

t(42,3)  -  C (41 , 2)  *  p lq 2  *  q4(3)  *  q3(2,3,3) 

P( 3, 2 )  -  p( 3, 2)  +  c ( 42 , 3 ) 
pni  tlt2c3  0  Table  6  Formula  3 

C(43,3)  =  t ( 41 , 2)  *  p2q 1  *  q4(j)  *  q3(2,3,3) 


p( 3,2)  -  p( 3,2  )  +  t ( “5, 3  ) 


t3  clt2t3  0 


Table  6 


Formula  4 


c 

t(44,J)  =  t ( 41 , 2 )  *  p 1P2  *  q4( 3 )  *  q3(2,3,3) 
c 

PC  3, 2)  =  p ( 3 , 2)  +  L ( 44 , 3 ) 

c  t2  tlc2  0  Table  6  Formula  1 

c 

t(45, 3)  »  t(45, 2)  *  qlq2  *  q4(3)  *  q3(l,3,2)  *  q3(2,3,3) 
c 

P( 3,2)  -  p( 3, 2)  +  t( 45, 3) 

c  t2t3  tlt2  0  Table  6  Formula  2 

c 

t ( 46, 3)  -  t(45, 2)  *  plq2  *  q4(3)  *  q3(l,3,2)  *  q3(2,3,3) 
c 

P( 3 , 2 )  -  p( 3, 2)  +  t ( 46 , 3 ) 

c  t2  cic2t3  U  Table  6  Formula  3 

c 

t(47, 3)  =  t ( 45 , 2)  *  p2q 1  *  q4(3)  *  q3(i,3,2)  *  q3(2,3,J) 
c 

P( 3, 2  )  -  p ( 3, 2 )  +  t ( 47 , 3) 

c  t2t3  tlc2t3  0  TaDle  o  Formula  4 

c 

t(48,3)  -  t ( 45 , 2 )  *  plp2  *  q4( 3)  *  q3(i,3,2)  *  q3(2,3,3) 
c 

PC 3, 2 )  *  p(3, 2)  +  t(48, 3) 

c  tl  tl  0  Table  6  Formula  1 

c 

t(49,3)  =  £ ( 49 , 2)  *  q lq 2  *  q4(3)  *  qJ(i,J,l)  *  q3( 2, 3,1) 
c 

P( 3 , 2 )  -  p ( 3 , 2 )  +  c ( 49 , 3) 

c  tlt3  tl  0  Table  6  Formula  2 

c 

t(50,3)  =  t ( 49, 2)  *  p lq 2  *  q4(3)  *  q3(l,3,l)  *  q3(2,3,l) 
c 

P( 3,2)  -  p( 3, 2)  +  t ( 50, 3) 


tl  tlt3  0  Table  6  Formula  3 

C(51,3)  -  t ( 49 , 2 )  *  P2q 1  *  q4(3)  *  q3(l,3,l)  *  q3(2,3,l) 
P( 3, 2)  -  PC  3, 2)  +  t( 51 ,3) 
clt3  tlt3  0  Table  6  Formula  4 

t ( 52 , 3 )  »  t(49,2)  *  p lp2  *  q 4 ( 3 )  *  q3(l,3,l)  *  q 3 ( 2 , 3 , 1 ) 
P( 3, 2)  -  p ( 3 , 2  )  -t-  t( 52, 3; 


r _ i  ~  \ 


[•I#*1! 


t(53,j)  =  t ( 53, 2)  *  qlq2  *  q4(3)  *  q3(l,3,3;  *  q3(2,3,l) 
P( 3, 2 )  -  p ( 3, 2 )  +  t ( 33,3) 


c  tlt2r3  cl  0  Table  b  Formula  2 

c 

t(  54, 3)  =■  C  (  53 , 2 )  *  plq2  *  q4(  3)  *  q3(l,3,3)  *  q3(2,3,l) 
c 

P(3, 2)  -  p( 3, 2)  +  t ( 54 , 3 ) 

c  clc2  tit 3  0  Table  6  Formula  3 

c 

t ( 55 , 3 )  =  t ( 53 , 2 )  *  p2q 1  *  q4(3)  *  q3(l,3,3)  *  q3(2,3,l) 
c 

P( 3,2j  =  p( 3, 2)  +  c( 55,3) 

—  — 

c  tlt2t3  tit3  0  Taole  o  Formula  4 

c 

t(5o, J)  =  t ( 53, 2)  *  p lp2  *  q 4 ( 3 )  *  q3(l,3,3)  *  q3(2,3,l) 
c 

P(j,2)  -  p (3,2)  +  t( 56, 3) 

c------------------------------ 

c  tl  tit 2  0  Table  6  Formula  1 

c 

t  (.  57 , 3)  =*  t  (  57, 2)  *  q  lq  2  *  q4(3)  *  q3(l,3,l)  *  q3(2,3,3) 
c 

P( 3, 2)  -  p( 3, 2)  +  t ( 57 , 3) 

c  tlt3  tlt2  0  Table  6  Formula  2 

c 

t  (  58 , 3)  =  t(57,2)  *  p  lq  2  *  q4(.3)  *  q3(l,3,l)  *  q3(2, 3,3) 
c 

P( 3, 2 )  =  p( 3, 2)  +  t( 58, 3) 

c  - 

c  tl  tlt2t3  0  Table  6  Formula  3 

c 

t( 59,3)  =  t ( 57 , 2)  *  p2q 1  *  q4(3)  *  q3(l,3,l)  *  q3(2,3,3) 
c 

p( 3, 2  )  =  p ( 3 , 2  )  +  t ( 59, 3) 

c  tlt3  tlt2t3  0  Table  6  Formula  4 

c 

t ( 60 , 3 )  »  t  C  57 ,2)  *  p lp2  *  q4(3)  *  q3(l,3,l)  *  q3(2,3,3) 
c 

P( 3, 2 )  -  p(3,2)  +  t ( 60, 3) 

c  tlt2  tlt2  0  Table  6  Formula  1 

c 

t(6i,3)  =  t ( 6 1 , 2 )  *  q lq 2  *  q4(3)  *  q3( 1 , 3 , 3 )  *  q3(2,3,3) 
c 

P( 3 , 2 )  =  p ( 3,2)  +  t ( 61 , 3 ) 

c  tlt2t3  tlt2  U  Table  6  Formula  2 

c 


77 


t(62,3)  =  t(61,2)  *  plq2  *  q4(3)  *  q3(l,3,3)  *  q3(2,3,3) 
P( 3, 2)  =  p( 3 , 2)  +  t  ( 62,  3) 
tlc2  clc2t3  0  Table  6  Formula  3 
t(63,3)  -  £(61.2)  *  p2q 1  *  q4(3)  *  q3(l,3,3)  *  q3(2,3,3) 
p( 3, 2)  -  p(3,2)  +  c(63,3) 


tlt2t3  t It2t3  0  Table  6 


Formula  4 


t(64,3)  -  t(61,2)  *  plp2  *  q4( 3)  *  q3(l,3,3)  *  q3(2,3,3) 
p(3,2)  -  p(3, 2)  +  t ( 64, 3) 


Table  b 


Formula  5 


t(65, 3)  = 


q4( 3)  *  ( 

+  t(i,; 
+  t(5,; 
+  t(9,: 
+  t( 13, 
+  t(17, 
+  c ( 21 . 
+  t  (25 
+  t(29, 
+  t(33 
+  t  ( 37 
+  t  ( 4 1 
+  C  ( 45 
+  t  ( 49 
+  t  (  53 
+  t(  57 
+  t  ( 6 1 


c ( b5 , 2) 


,2)  *  (  1.  - 

,2)  *  (  1.  - 

,2)  *  (  1.  - 

,2)  *  (  1.  - 

,2)  *  (  1.  - 

,2)  *  (  1.  - 

,2)  *  (  1.  - 

,2)  *11.- 
,2)  *  (  1.  - 

,2)  *  (  1.  - 

,2)  *  (  1.  - 

.2)  *  (  1.  - 

.2)  *  (  1.  - 


*  q3( 1 , 3 , 2) 

*  q3(  2 , 3 , 2 ) 

*  q3( 1 , 3 , 2 ) 

*  q3(l,3,l) 

*  q3(l,3,3) 

*  q3( 1,3,1) 

*  q3( 1 , 3, 3) 

*  q  3 ( 2,3,1) 

*  q 3 C  1,3,2) 

*  q3(2,3,3) 

*  q3( 1 , 3 , 2 ) 

*  q 3 ( 1,3,1) 

*  q3( 1 , 3 , 3) 

*  q3(  1 , 3 , 1 ) 

*  q3(l,3,3) 


*  q3(2,3,2)  ) 

) 

) 

*  q3(2,3,2)  ) 

*  q3( 2, 3 , 2)  ) 

) 

*  q3( 2 , 3 , 1 )  ) 

) 

*  q3( 2, 3 , 3 )  ) 

*  q3(2,3,l)  ) 

*  q3(2,3,l)  ) 

*  qj(2,3,3)  ) 

*  q3(2,3,3)  )  ) 


P( 3, 4 )  =  p(3,4)  +  t ( 65, 3 ) 
P(3,8)  =  p(3,8)  +  t(65,3) 
P(3,9)  «  p( 3, 9 )  +  t ( 65 , 3) 


phi  phi  1 


Table  6  Formula  6 


C ( 1,6)  -  q lq2  *  (  t( 1.2)  *  p4(3)  +  t(l,5)  ) 

P( 3, 3)  -  p( 3, 3)  +  t ( 1 , 6) 
p( 3 , 7 )  -  p( 3 , 7 )  +  t ( 1,6) 
p ( 3,9)  -  p( 3, 9)  +  t ( 1 , 6) 


t3  phi  1 


Table  o  Formula  7 


t(2,b)  =  p lq 2  *  (  t(l,2)  *  p4(3)  +  c ( 1,5)  ) 


p(  3, 5)  *  p(  3, 5 )  -r  c ( 2 , 6) 
p( 3 , 7 )  =  p( 3,7)  +  t(2,6) 

P( 3, 9)  -  p( 3, 9)  +  t ( 2, 6 ) 
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p(3,9)  +  t (9,6) 


c 

c 

c 

c 


c 

c 

c 

c 


c 

c 

c 

c 


c 

c 

c 


c 


c 

c 

c 


c 


c 

c 

c 


c 


t3  c2  1  Table  b  Formula  7 

t ( 10, b)  -  plq2  *  q3U,J,2)  *  (  t(9,2)  *  p4(3)  +  t(9,5)  ) 

p(3,5)  -  p(3, 5)  +  t( 10,6) 

P(3,7)  -  p<3,7)  +  t( 10,6) 

P( 3, 9)  -  p( 3, 9)  +  t( 10, 6) 


phi  t2t3  1  Table  6  Formula  8 

t( 11,6)  -  p2q 1  *  q3(2,3,2)  *  (  t(9,2)  *  p4(3)  +  t(9,5)  ) 

p(  3, 5)  =  p( 3, 5)  +  t ( 1 1 , 6) 
p( 3 , 7 )  -  p( 3 , 7 )  +  t ( 1 1 , 6) 
p(  3, 9)  -  p (  3 , 9 )  -r  t ( 1 1 , 6) 


c3  c2t3  1  Table  6  Formula  9 

t(12,b)  -  plpZ  *  q3(2,3,2)  *  (  t(9,2)  *  P4(3)  +  t(9,5)  ) 

p(  3,5)  =■  p(  3, 5)  +  c( 12, 6) 
p( 3 , 7)  =  p( 3 , 7 )  +  t( 12,6) 
p( 3, 9)  -  p( 3, 9)  +  c( 12,6) 


t2  t2  1  Table  6  Formula  6 

C ( 13,6)  -  q lq2  *  q3(x,3,2)  *  q3(2,3,2) 

*  (  t( 13,2)  *  p4(3)  +  c ( 13,5)  ) 

p(  3, 5)  =*  p(  3,  5)  +  c(  13,6) 
p( 3 , 7 )  -  p( 3 , 7 )  +  t ( 13,6) 
p( 3 , 9)  -  p( 3, 9)  +  t ( 13 , 6) 


t2c3  c2  1  Table  6  Formula  7 

ECU,  6)  -  p  lq  2  *  q  3(  1 , 3 , 2 )  *  q3(2,3,2) 

*  (  t(  13,2)  *  p4( 3)  +  t ( 13,5)  ) 

P( 3, 5)  -  p( 3, 5)  +  t ( 14, 6) 
p( 3, 7)  -  p ( 3,7)  +  t ( 14,6) 
p( 3, 9)  -  p( 3, 9)  +  t( 14, b) 


t2  t2t3  1  Table  6  Formula  8 

1(15,5)  -  P2q 1  *  q3( 1,3,2)  *  q3(2,3,2) 

*  (  t( 13,2)  *  p4( 3 )  +  t( 13,5)  ) 

p( 3, 5)  -  p( 3, 5)  +  c ( 13,6) 
p( 3, 7)  -  p( 3 , 7 )  +  t ( 15,6) 
p( 3,9)  -  pC  3, 9)  +  t( 15, b; 


c 

c 

c 


c2t  3  1 2t  5  1 


Table  6 


Formula  9 


t(22,6)  -  plq2  *  q 1 , 3, 3) 


c 


*  (  t ( 21 , 2)  *  p4(  3)  +  t(21 ,5)  ) 


p(i,5)  =  p  ( 3 , 5 )  +  t ( 22, b) 

P( 3, 7)  «  p( 3 , 7 )  +  t ( 22 , b ) 
p( 3, 9)  =  p<.3,9)  +  t ( 22, b) 

c  tlt2  t3  1  Table  6  Formula  8 

c 

t(23,6)  -  P2q 1  *  q3( 1 ,3,3) 

*  *  (  t(21,2)  *  p4(  3)  +  t(21,5)  2 
c 

P( 3, 5)  -  p( 3, 5)  +  t ( 23, 6) 

P( 3 , 7)  -  p(3,7)  +  t(23,6) 
p( 3, 9)  -  p( 3,9)  +  c(23,6) 

c  tlc2t3  t3  1  Table  6  Formula  9 

c 

t(24,o)  =  plp2  *  q3( 1 , 3 , 3 ) 

*  *  (  t (21 ,2)  *  p4( 3)  +  c (21 ,5)  ) 
c 

p( 3, 5)  -  p( 3, 5)  +  t ( 24 , 6) 
p( 3, 7)  -  p( 3, 7)  +  t(24,6) 

P( 3, 9)  -  p(3,9)  +  t ( 24 , 6 ) 


c  cl  c2  1  Table  6  Formula  6 

c 

t(25,6)  -  qlq2  *  q3( 1 , 3 , 1 )  *  q3(2,3,2) 

*  *  (  t ( 25 ,2)  *  P4(3)  +  t ( 25 , 5)  ) 
c 

p(  3, 5)  =■  p(  3, 5 )  +  t(25,6) 
p(3,7)  -  p( 3, 7)  +  t(25,6) 

P( 3, 9)  «  p(  3, 9)  -r  t(25,6) 

c  tit 3  t2  1  Taole  b  Formula  7 

c 

t(2b, 6)  =  p lq 2  *  q3( 1 ,3 , 1 )  *  q3(2,3,2) 

*  *  (  t(25,2)  *  P4( 3 )  +  t (25 ,5)  ) 
c 

p( 3, 5)  -  p( 3, 5 )  +  t ( 2b, 6) 
p(3,7)  -  p( 3 , 7 )  +  t(2b,b) 

P( 3, 9)  -  p( 3,9)  +  t ( 26, 6) 

c  tl  t2t3  1  Table  6  Formula  8 

c 

t (27, 6)  -  p2q 1  *  q3( 1,3,1)  *  q3(2,3,2) 

*  *  (  t(25,2)  *  p4( 3 )  +  t(25,5)  ) 
c 

P( 3, 5)  -  p( 3, 5)  +  t ( 27 , 6) 
p( 3 , 7)  -  p( 3 , 7 )  +  t ( 27 , 6 ) 

P( 3, 9)  -  p ( 3 , 9 )  +  t ( 27 , 6 ) 


c  tlt3  t2t3  1  TaDle  6  Formula  9 

c 

t ( 28 , 6)  =  plP2  *  q 3(  1 , 3, 1 )  *  q?(2,3,2) 

*  *  (  t(25,2)  *  P4(3)  +  c ( 25 , 5)  ) 
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At 


tit  2  t2  i 


Table  6 


Formula  6 


t(29,6)  -  q lq2  ■*  q3(l,3,3)  *  q3(2,3,2) 

*  (  t(29,2)  *  p4(3)  +  t(29,5)  ) 

P(3,5)  *  p( 3, 5)  +  t(29,6) 

P(3,7)  -  p(3,7)  +  t(29,6) 
p(3,9)  -  p( 3,9)  +  t (29 , 6) 


tlt2t3  t2  1  Table  6  Formula  7 

t(30,6)  *  Plq2  *  q3(l,3,3)  *  q3(2,3,2) 

*  (  t(29,2)  *  p4(3)  +  t(29,5)  ) 

p( 3, 5)  =*  p( 3, 5)  +  t ( 30, 6) 
p( 3, 7 )  =  p( 3, 7)  +  t(30,6) 

P(3,9)  =  p(3,9)  +  t ( 30, 6) 


tlt2  t2t3  1  Table  6  Formula  8 

t(31 ,6)  -  P2q 1  *  q3( 1 , 3 , 3)  *  q3(2,3,2) 

*  (  t(29,2)  *  p4(3)  +  t(29,5)  ) 

P( 3, 5)  -  p( 3, 5)  +  t ( 31 , 6) 
p(3,7)  -  p(3,7)  +  t(31 ,6) 

P( 3, 9)  -  p( 3, 9)  +  t(31,6) 


tlt2t3  t2t3  1  Table  6  Formula  9 

t(32,6)  =*  p 1P2  *  q3( 1,3,3)  *  q3(2,3,2) 

*  (  t(29,2)  *  P4( 3)  +  t(29,5)  ) 

P( 3, 5)  -  p( 3,5)  +  t(32,6) 

P( 3, 7)  -  p( 3, 7)  +  t(32,6) 

P( 3, 9)  -  p( 3, 9 )  +  t(32,6) 


phi  tl  1  Table  6  Formula  b 

t(33,6)  -  qlq2  *  q3(2,3,l) 

*  (  t(33, 2)  *  P4(3)  +  t(33,5)  ) 

P( 3, 5)  -  p( 3, 5)  +  t ( 33, 6) 

P( 3, 7)  -  p( 3 , 7 )  +  t ( 33 , b) 
p( 3,9)  -  p(3,9)  +  t ( 33 , 6 ) 


t3  tl  1  Table  b  Formula  7 

t(34,b)  *  plq2  *  q3( 2, 3 , 1 ) 

*  (  tl 33, 2)  *  p4(3)  +  t ( 33 , 5 )  ) 


phi  tlc3  1 


Table  6 


Formula  8 


c 


c 


c 

c 

c 


c 


c 

c 

c 


c 


c 

c 

c 


c 


c 

c 

c 


c 


c 

c 

c 


c 


t(35, 6)  -  P2ql  *  q3(2, 3, 1 ) 

*  *  (  t(33,2)  *  p4(  3)  +  t ( 33, 5)  ) 

P(3,5)  -  p( 3, 5)  +  t ( 35, 6) 

P(3,7)  -  p( 3, 7)  +  t( 35 ,6) 

P( 3, 9)  -  p( 3, 9)  +  t(35, 6) 


Table  6  Formula  9 
*  q3( 2 , 3, 1 ) 

*  (  c ( 33 , 2 )  *  p4(3)  +  t(33,5)  ) 

P( 3, 5)  =  p( 3, 5)  +  t ( 36 , 6) 

P( 3, 7)  =  p( 3, 7)  +  c(36,6) 

P( 3, 9)  =■  p( 3, 9)  +  t( 36, 6) 


t2  tl  1  Table  6  Formula  b 

t(37,6)  -  qlq2  *  q3(l,3,2)  *  q3(2,3,l) 

*  *  (  t(37 ,2)  *  P4(3)  +  C ( 37 , 5)  ) 

p(3,5)  *  p( 3, 5)  +  t ( 37 , 6) 

P( 3 , 7)  -  p( 3 , 7)  +  t( 37 ,6) 
p(3,9)  «  p( 3,9)  +  t ( 37 , 6) 


t2t3  tl  1  Table  6  Formula  7 

t(38,6)  *  plq2  *  q3(l,3,2)  *  q3(2,3,l) 

*  *  (  t(37,2)  *  P4( 3)  +  t(37 ,5)  ) 

p( 3 , 5)  -  p(3, 5)  +  t ( 38, 6) 
p( 3 , 7)  -  p( 3, 7)  +  t(38,6) 
p( 3, 9)  -  p( 3, 9)  +  t(38,6) 


t2  tlt3  1  Table  6  Formula  8 

t(39,6)  -  P2q 1  *  q3( 1,3,2)  *  q3(2,3,l) 

*  *  (  t(37 ,2)  *  P4( 3)  +  t(37 ,5)  ) 

P( 3, 5)  -  p( 3, 5)  +  t ( 39 , 6) 

p( 3 , 7)  -  p( 3,7)  +  t ( 39, 6) 
p( 3, 9)  -  p(3,9)  +  t(39 ,6) 


t2t3  tlt3  1  TaDle  6  Formula  9 

t(40,6)  -  p ip2  *  q3(  1, 3,2)  *  q 3 <. 2 , 3 , 1) 

*  (  c( 37,2)  *  P4(3)  +  t( 37,5)  ) 


t3  t lt3  1 
t(36,b)  *  plp2 


P( 3, 5)  -  p( 3, 5)  +  c(40,6) 


pni  tlt2  1  Table  6  Formula  6 

t(41,6)  -  qlq2  *  q3(2,3,3) 

*  (  t(41 , 2)  *  P4( 3)  +  t(41,5)  ) 

p(3,5)  -  p( 3, 5)  +  t(41,6) 
p(3,7)  -  p( 3,7)  +  t(41,6) 

P( 3, 9)  -  p( 3, 9)  +  t(4i, 6) 


1  Table  6  Formula  7 

plq2  *  q3(2, 3, 3) 

*  (  c(41 ,2)  *  P4( 3)  +  t (41 , 5)  ) 

p( 3, 5)  =  p(3,5)  +  c(42,b) 
p(3,7)  *  p( 3, 7 )  +  t(42,6) 
p(3,9)  =  p(3,9)  +  t(42, 6) 


phi  clc2t3  1  Table  6  Formula  8 

c(43,6)  -  P2q 1  *  q3( 2, 3 , 3) 

*  (  t (41 , 2)  *  P4(3)  +  c(41,5)  ) 

P( 3 , 5 )  -  p(3,5)  +  t (43, 6) 
p(3,7)  ■  p  ( 3 , 7 )  +  c(43,6) 

P( 3,9)  -  p( 3, 9)  +  C ( 43 , 6 ) 


t3  t  it 2 
t ( 42, 6)  - 


t3  tlt2c3  i  Table  b  Formula  9 

t(4<*,6)  =  p ip2  *  q3(2,3,3) 

*  (  c(41 ,2)  *  p4 ( 3 )  +  t (41 , 5)  ) 

p(  3,5)  -  p( 3, 5)  +  t(44,6) 

P(  3,7)  -  p( 3 , 7)  +  t(44,6) 

P( 3, 9)  =  p ( 3 , 9 )  +  C(44,b) 


t2  tlt2  1  Table  b  Formula  b 

t(45,6)  -  qlq2  *  q3(l,3,2)  *  q3(2,3,3) 

*  (  t(45,2)  *  P4( 3)  +  t(45,5)  ) 

P( 3, 5)  -  p( 3, 5)  +  t(45,b) 
p( 3 , 7)  -  p( 3,7)  +  t(45,6) 

P(3,9)  -  p(3,9)  +  t ( 45, b ) 


t2t3  tlt2  1  Table  b  Formula  7 

C(4b,b)  -  plq2  *  q3( 1 , 3, 2)  *  q3(2,3,3) 

*  (  t  ( 45 , 2 )  *  p<*(3)  +  t  ( 45 , 5)  ) 


P ( 3 , 5 )  ■  p( 3, 5)  +  t(4b, 6) 
p( 3, 7)  -  p( 3 , 7 )  +  t(4b,b) 


1 


f 

r 

f 

I 


! 


k 

L 


\ 


c  c2  tlc2c3  i  Table  b  Formula  b 

c 

t(47,o)  -  P2q 1  *  q3( 1 , 3 , 2)  *  q3(2,3,3) 

*  *  C  t(45 ,2)  *  p4(3)  +  t(45,5)  •) 

c 

p(3,5)  -  p( 3, 5)  +  t (47 , 6) 

P( 3 , 7)  -  p(3,7)  +  t( 47 ,6) 

P( 3, 9)  -  p( 3,9)  4-  t(47,6) 


c  t2t3  tlt2t3  1  Table  6  Formula  9 
c 

t ( 4b , 6 )  -  p Ip2  *  q3(l,3,2)  *  q3( 2, 3, 3) 

*  *  (  t(45,2)  *  P4(3)  +  t(45,5)  ) 

c 

P( 3, 5 )  =  p( 3, 5)  +  c(4b, 6) 
p( 3 , 7)  *  p( 3 , 7)  +  t(4b,6) 
p( 3, 9)  -  p( 3, 9)  +  t ( 4b , 6) 


c  tl  tl  1  Table  6  Formula  b 


t ( 49 , 6 )  -  q lq2  *  q3(l,3,l)  *  q3(2,3,i) 

*  *  (  t(49,2)  *  p4(3)  +  t(49,5)  ) 
c 

P( 3, 5)  -  p(3, 5)  +  t ( 49, 6) 

P( 3, 7)  «  p( 3, 7)  +  t(49,6) 

P( 3 , 9)  -  p(3,9)  +  C ( 49, 6) 
c----------------------- 

c  tlc3  tl  1  Table  6  Formula  7 

c 

c( 50, 6)  -  p lq 2  *  q 3 ( 1 , 3 , 1 )  *  q3(2,3,l) 

*  *  (  t(49,2)  *  p4(3)  +  t(49 ,5)  ) 
c 

P( 3, 5)  -  p(3,5)  +  t(50,6) 

P( 3,7)  =  p( 3 , 7)  +  c( 50,6) 
p( 3,9)  *  p( 3, 9)  +  t ( 50,6) 

c  cl  clt3  1  Table  6  Formula  b 

c 

t ( 51 , b )  -  P2q 1  *  q 3( 1 , 3 , 1)  *  q3(2,3,l) 

*  *  (  t(49,2)  *  P4( 3)  +  t(49,5)  ) 
c 

P( 3, 5)  -  p( 3, 5)  +  t ( 51 , 6) 

p( 3, 7 )  -  p( 3 , 7)  +  t( 51 ,6) 

P( 3, 9)  -  p(3,9)  +  t ( 51 , 6) 

c----------------------- 

c  tlt3  Clt3  1  Table  6  Formula  9 


t  (  52 , 6 )  -  p  1  p  2  *  q3U,3,l)  *  q3(2,3,l) 

*  (  t ( 49 , 2 )  *  p4( 3)  +  t(49,5)  ) 

p( 3, 5)  -  p( 2, 5)  +  c( 52,6) 
p( 3, 7)  -  p( 3 , 7 )  +  c( 52,6) 
p( 3, 9)  -  p( 3, 9)  +  c ( 52, 6) 


86 


t  lt2  tl  1 


Table  6 


Formula  6 


t(5J,6)  =  qlq2  *  q3( 1 , 3 , 3 )  *  q3(2,3,l) 

*  *  (  t( 53,2)  *  p4(3)  +  C ( 53 , 5)  ) 

P( 3, 5)  -  p(3,5)  +  t( 53, 6) 

P( 3, 7)  -  p( 3, 7)  +  t(53,6) 

P(3,9)  -  p( 3, 9)  +  t(53,6) 


tlt2t3  tl  1  Table  6  Formula  7  _ 

t(54,6)  -  Plq2  *  q3( 1 ,3, 3)  *  q3(2,3,l) 

*  (  t( 53,2)  *  P4(3)  +  t(53,5)  ) 

p( 3, 5)  -  p( 3,5)  +  t ( 54 , 6 ) 

P( 3, 7)  -  p( 3 , 7)  +  t(54,6) 

P( 3,  9)  *  p(  3,  9)  -r  t  ( 54, 6) 


tlt2  tlt3  1  Table  6  Formula  8 

t(55,t>)  =  P2q  1  *  q3(  1,3,3)  *  q3(2,3,l) 

*  *  <  t(53,2)  *  P4( 3)  +  t(53,5)  ) 

p( 3,5)  *  p( 3, 5)  +  t ( 55, 6) 
p( 3 , 7)  -  p( 3 , 7 )  +  t(55,6) 
p( 3, 9)  *  p( 3,9)  +  t(55,6) 


tlt2t3  tit3  1  Table  6  Formula  9 

t(5b,6)  =  p lp2  *  q3( 1,3,3)  *  q3(2,3,l) 

*  (  t(53,2)  *  P4(3)  +  t ( 53 , 5 )  ) 

p( 3, 5)  *  p(3,5)  +  t ( 56 , 6 ) 

P( 3, 7)  -  p( 3, 7)  +  t( 56, 6) 

P( 3, 9)  -  p( 3, 9)  +  t(56,6) 


tl  tlt2  1  Table  6  Formula  b 

t( 57 , 6)  -  qlq2  *  q3(i,3,l)  *  q3(2,3,3) 

*  (  t ( 57 ,2)  *  P4(3)  +  t ( 5 7 ,5)  ) 

p( 3, 5)  -  p( 3, 5)  +  t ( 57 , 6) 
p( 3, 7)  -  p( 3, 7)  +  t( 57,6) 
p( 3, 9)  -  p( 3, 9)  +  t ( 57 , 6) 


tlt3  tit 2  1  Table  6  Formula  7 

tC5B,6)  -  plq2  *  q3(l,3,l)  *  q3(2,3,3) 

*  (  t( 57 ,2)  *  p4( 3)  +  t( 57,5)  ) 


p( 3, 5)  -  p( 3, 5)  +  t( 58,6) 
p( 3, 7)  -  p(3, 7)  +  t ( 58 , 6) 
p( 3, 9)  -  p( 3, 9)  +  t ( 58, 6 J 


tl  tlt2t3  1  Table  6  Formula  6 

t( 59 , 6)  =  p2q 1  *  q3( 1 , 3, 1 )  *  q3(2,3,3) 

*  (  t(57,2)  *  p4(3)  +  t(57 ,5)  ) 


p(3,5)  -  p(3, 5)  +  c(59 ,6) 
P(3,7)  -  p(3,7)  +  t(59,6) 
p( 3, 9)  -  p( 3, 9)  +  t(59,6) 


tlc3  tlt2t3  1  Table  6  Formula  9 


c(60,6)  *  plP2  *  q3( 1 , 3 , 1 )  *  q3(2, 3, 3) 

*  (  t( 57 ,2)  *  p4(3)  +  t ( 57 ,5)  ) 

P( 3, 5)  -  p( 3, 5)  +  t(60,6) 
p(3,7)  -  p(3,7)  +  t( 60,6) 
p( 3 , 9)  *  p  ( 3 , 9 )  +  t(60,6) 


Clt2  tic2  1  Table  6  Formula  6 


t ( bi , 6)  =  qlq2  *  q3(l,3,3)  *  q3(2,3,3) 

*  (  t(61 ,2)  *  p4(3)  +  t(61 ,5)  ) 

p(3, 5)  *  p(3,5)  +  t ( 61 , 6) 
p( 3 , 7)  »  p( 3, 7)  +  c(61 ,6) 

P(3,9)  -  p( 3, 9)  +  t(61 ,6) 


tlt2c3  tlt2  1  Table  6  Formula  7 

t ( 62 , b )  =  plq2  *  q3(l,3,3)  *  q3(2,3,3) 

*  *  (  c(61,2)  *  P4(3)  +  t ( 61 , 5)  ) 

P( 3, 5)  -  p( 3, 5)  +  t(62,6) 
p(3,7)  -  p( 3, 7)  +  t(62,6) 

P( 3, 9)  -  p( 3, 9)  +  t ( 62 , 6) 


tlt2  tlt2t3  1  Table  6  Formula  8 

t ( 63 , 6)  *  P2q 1  *  q3( 1 , 3, 3)  *  q3(2,3,3) 

*  (  c(61,2)  *  P4(3)  +  t(61,5)  ) 

p( 3, 5)  -  p( 3, 5)  +  t ( 63, 6) 
p(3,7)  -  p( 3 , 7)  +  t ( 63, 6) 
p( 3,9)  -  p( 3,9)  +  t ( 63 , 6) 


tlt2t3  tlt2t3  1  Table  6  Formula  9 

t ( 64 , 6 )  -  plp2  *  q3( 1,3,3)  *  q3(2,3,3) 

*  (  t( 61 ,2)  *  P4(3)  +  t(61 ,5)  ) 


p( 3 , 5 )  «  p( 3, 5)  +  t ( 64 , 6 ) 
p( 3, 7)  -  p( 3, 7)  +  t(64,b) 
p(3,9)  -  p( 3, 9;  +  t ( 64, 6) 


t(65,b)  *  t(65,2)  *  p 4(3)  +  c(65,5) 

*  (  c ( 1,2)  *  p4(3)  +  t(l,5)  )  *  ( 

k  +  (  t(5,2)  *  P4(3)  i-  t(5,5)  )  *  ( 

*  +  (  t(9,2)  *  P4(3)  +  t(9,5)  )  *  ( 

*  +  (  c(l 3,2)  *  P4(3)  +  t(13,5)  ) 

*  *  (  1.  -  q5q6  * 

*  +  (  t ( 17 , 2)  *  P4(3)  +  t(  17 , 5)  )  * 

*  4-  (  t(21 ,2)  *  p4(3)  +  t(21 ,5)  )  * 

*  +  (  t ( 25 , 2)  *  p4(3)  +  t(25,5)  ) 

*  *  (  i.  -  q5q6  * 

*  +  (  t ( 29 , 2)  *  P4( 3)  +  t ( 29 , 5 )  ) 

k  *(1.  -q5q6* 

*  +  (  t(33,2)  *  P4( 3)  +  t ( 33, 5)  )  * 
t(65,6)  -  c(65,6)  +  (  t(37,2)  *  p4(3) 

*  *  (  1.  -  q5q6  * 

*  +  (  c ( 41 , 2)  *  P4( 3)  +  t(41 ,5)  )  * 

*  +  (  t(45, 2)  *  p4( 3;  +  t ( 45, 5 )  ) 

*  *(1.  -q5qb* 

*  +  (  t(49,2)  *  Ph(3)  +  t ( 49 , 5 )  ) 

*  *  (  1.  -  q5qo  * 

*  +  (  t ( 53, 2)  *  P4(3)  +  t(53,5)  > 

*  *  (  1.  -  q5q6  * 

*  +  (  t ( 57 , 2)  *  P4( 3)  +  t ( 57 , 5 )  ) 

*  *  (  1.  -  q5q6  * 

*  +  (  t(61 ,2)  *  P4( 3)  +  £(61,5)  ) 

*  *  (  1.  -  q5qb  * 


1.  -  q5q6  ) 
i .  -  q5qb  * 
1.  -  q5q6  * 

q3( 1,3,2)  * 
(  1.  -  q5q6 
(  1.  -  q5q6 


q3( 1 , 3 , 2 )  ) 
q3(2, 3,2)  ) 

q  3  ( 2 , 3 , 2 )  ) 

*  q3( 1,3,1)  ) 

*  q  3 ( 1,3,3)  ) 


q 3 ( 1,3,1)  *  q3(2,3i2_)  ) 


q3( 1,3,3)  * 
(  1 .  -  q5q6 
+  c  C  37 , 5 ) 
q3( 1,3,2)  * 
(  1.  -  q5q6 


q3( 2,3,2)  ) 

*  q3( 2 , 3 , 1 )  ) 

) 

q 3(  2 , 3 , 1 )  ) 

*  q  3(  2 , 3 , 3 )  ) 


q 3 C  1,3,2)  *  q 3( 2 , 2 , 3)  ) 
q  3( 2 , 3 , 1 )  *  q3(2,3, 1)  ) 
q 3 ( 1,3,3)  *  q3( 2 , 3 , 1 )  ) 
q3( 1 , 3 , 1 )  *  q3(2 , 3 , 3)  ) 
q 3( 1,3,3)  *  q3( 2 , 3 , 3 )  ) 


pC 3, 6)  *  p( 3, 6)  +  t (b5, 6) 
P ( 3 , 7 )  «  p( 3 , 7 )  +  t(65,b) 
P( 3,8)  »  p( 3,8)  +  t(b5, 6) 
p(3,9)  =*  p(  3 , 9 )  +  t(65,6) 


ptot  -  p( 3, 1 )  +  p( 3 , 2 )  +  p( 3,3)  +  p( 3 ,4 )  +  p(3,5)  +  p(3,6) 
if (  abs(  1.  -  ptot  )  .gt.  .000001  )write  (*,100)  ptot 
100  format('  sum  of  INDEPENDENT  probabilities  for  3  bursts  is  ' ,fl0.7) 
return 
end 
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